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ABSTRACT 

o ■ 

in: 

c/3 , We present a simplified chemical and thermal model designed to allow com- 

I— I, putationally efficient study of the thermal evolution of metal-poor gas within 

^ ! large numerical simulations. Our main simplification is the neglect of the molec- 

^ ! ular chemistry of the heavy elements. The only molecular chemistry retained 

00 '. within the model is the formation and destruction of molecular hydrogen. De- 

^ i spite this major simplification, the model allows for accurate treatment of the 

', thermal evolution of the gas within a large volume of parameter space. It is valid 

O I for temperatures 50 < T < 10000 K and metallicities < Z < 0.1 Zq. In gas 

Q I with a metallicity Z = 0.1 Z©, and in the absence of an incident ultraviolet radi- 

^ I ation field, it is valid for hydrogen number densities nn ^ 500/tchar cm~^, where 

I tchar is the sizc in Myr of the characteristic physical timescale of interest in the 

^ I problem. If Z -C 0.1 Zq, or if a strong ultraviolet radiation field is present, then 

■ " " ' the model remains accurate up to significantly higher densities. We also discuss 

some possible applications of this model. 

Subject headings: astrochemistry — molecular processes — ISM: molecules - 
galaxies: formation - cosmology: theory 



1. Introduction 



It has long been known that cooling by molecul ar hydrogen, H2, plays a major role in 
regulating star formation in primordial protogalaxies (ISaslaw &: Zipoylll967l : [Peebles &: Dicke 



- 2 - 



19681 : iMatsuda et al.lll969l : iLepp fc Shul]||l983l : IXegmark et al.lll997l : lAbel et al.ll2002h . The 
importance of H2 stems from the fact that in most circumstances it is the dominant coolant in 
primordial gas at T < 10^ K. As a result, the cheraistry of in primord i al gas has attracted 



considerable study JPakarno fc LeDD lll987l:lBlackl ll99ll:lAbel etaDll997l : [Gam fc PallalE998 



Standi. Lepp. fc Dalgarnol 19981 : Lepp. Stancil. fc Dalgarnol 2OO2I ) and it is now generally 



accepted that only a small chemical network, of maybe 20-30 reactions, is required to model 
H2 chemistry over a very wide range of conditions in these systems. 

The introduction of metals into the gas, as will occur following enrichment of their 
surroundings by the first generation of supernovae, complicates matters enormously. Many 
other atomic and molecular coolants become available and the associated chemistry is highly 
complex: for instance, a reasonably complex model of purely gas-phase chemistry can eas- 
ily st retch to ~ 400 reactants and almost 4000 reactions (e.g. iLe Teuff. Millar fc Markwick 
2OOOI ) . This would not matter were it not for two important points. First, many astrophysi- 
cists believe that metal enrichment above a certain level - the so-called 'critical metallicity' 
- leads to a significant change in the stellar initial mass function (IMF), from an IMF domi- 
nated by massive stars to one that looks far more like the familiar Salpeter IMF of local star 
formation (IBromm et al.ll200ll : ISchneider et al.ll2002l ). Second, testing this idea numerically 
using three-dimensional hydrodynamical simulations requires us to model the chemistry of 
the gas, but the highly detailed chemical models mentioned above are impractical to use 
in high-resolution numerical simulations, owing to their high computational cost. This is a 
consequence of the fact that chemical rate equations are frequently stiff and so for reasons 
of stability must be solved implicitly, with a computational cost that scales as the cube of 
the number of chemical species involved. 

It is therefore important to look for ways to simplify the chemistry without unduly 
compromising the accuracy of the resulting model. Major simplifications can be made if we 
make the reasonable assumption that the main coolants in low metallicity, high redshift gas 
will be similar to those in the local interstellar medium - which should be true provided that 
the abundance ratios of the various metals in low metallicity gas are not too unusual - and 
also if we restrict the range of physical conditions in which we are interested. 

In this paper, we present a chemical model designed to model the chemistry of the 
major coolants in cool (T < 10'^ K) low-metallicity gas at low gas densities. In the ab- 
sence of ultraviolet radiation, and in gas with Z = 0.1 Z0, our model is valid for atomic 
hydrogen number densities nn ^ 500/tchar cm"'^, where tchar is the size in Myr of the char- 
acteristic physical timescale of interest in the problem, corresponding to an overdensity of 
5 = 2.5xl0^t~^^^{l+z)~^ with respect to the cosmological mean background density. At lower 
metallicities, or if a moderately strong UV field is present, our model remains accurate up to 
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significantly higher densities. We do not treat the coohng or c hemistry of hot gas (?^^ 10 K) 



as this has aheady been treated in detail elsewhere (see e.g. [Sutherland fc Dopitalll993l ). 



The structure of this paper is as follows. In §2] we present our simplified chemical 
network and discuss the principles determining our choice of reactants ( §2.11) and reactions 
( §2.2p . as well as our treatment of grain surface chemistry ( §2.31) and photochemistry ( §2.4p . 
In ^ we discuss our treatment of the main thermal processes included in our model, with 
a particular emphasis on atomic fine structure cooling. We conclude in §1] with a brief 
discussion of possible applications of our model. 



2. Chemical model 



2.1. Choice of chemical species 



Our choice of which chemical species to include in our simplified network was guided 
by two main considerations. In order to properly model the thermal evolution of the gas, 
we must be able to accurately model the evolution of the chemical abundances of all of the 
major coolants. At the same time, in order to keep our chemical treatment computationally 
efficient, we do not want to include more species than are strictly necessary. In some cases, 
the decision on whether or not to include a species was obvious. For instance, neutral atomic 
hydrogen, H, is a major constituent of the gas and is also a major coolant at T ~ 10^ K 
and above. Molecular hydrogen, H2, also must be included as it has long been recognised 



to be the domiriant co olant in metal-free gas at 200 < T < 10^ K ( jSaslaw fc Zipoy 



1967 



Peebles fc Dickd 119681). and it is now clear that it also remains si gnificant in metal-poor 
gas (lOmukai et al.ll2005l : ISantoro fc Shulll l2006l : IJappsen et al.l 120071 ). A good case can also 
be made for the inclusion of H P, which dominates the cooling in zero-met allicity gas at 
T < 200 K Jpiower et allboooh. and which remains important at least up to metallicities 
Z ~ 10"^ Zq ( lOmukai et al.ll2005l ). As far as metals go, however, the choice is less obvious. 
Metal atoms and ions produce little in the way of resonance line cooling at T < 10^ K, but 
many can act as sources of fine structure line emission or emission from metastable transitions 
at these temperatures. Therefore, in order to determine which of the various species are 
important coolants in the region of parameter space that our model is designed to cover, we 
directly compared the cooling rates produced by each species, under the assumption that 
the relative abundances of the various elements were the same as in the local interstellar 
medium. 

The species we investigated in th i s corn parison were the same as those included in the 
shock models of iHoUenbach fc McKeel Jl989h : C, C+, CI, C1+, Fe, Fe+, N, N+, Ne+, Ni, Ni+, 
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O, O"*", S, S"*", Si and Si"*". This list includes all of the neutral or singly ionized species that 
have both non- negligible abundances and emission lines that are accessible at T < 10^ K. 
Doubly ionized species are unlikely to be abundant in gas of this temperature, and so are 
unlikely to contribute significantly to the total cooling rate. Atomic data for C, C+, O, Si 
and Si"*" was tak en from the sources lis t ed in Tables [5H6l For the remaining species, we used 
data taken from iHoUenbach fc McKed (119891). For the relative abu ndances of C, Fe, N, Ne, 
Ni, O, S and Si, w e used values taken from ISembach et al.l (120001 ). while for CI we used a 



value taken from lAsplund. Grevesse fc Sauvall (120061 ). 



We investigated three different scenarios for the ionization state of the gas, which we 
here refer to as the 'no ionization', 'moderate ionization' and 'high ionization' cases. In 
the no ionization case, we assumed that the gas was completely neutral, and so set the 
electron abundance to zero, along with the abundances of all of the ionized species. In the 
moderate ionization case, we assumed that the fractional ionization of hydrogen was such 
that nn+fn-n = 10~^, that species with ionization potentials lower than that of hydrogen - 
C, CI, Fe, Ni, S and Si - were fully ionized, and that all of the other species remained fully 
neutral. Finally, in the high ionization case, we assumed that ^h+Z'^.h = 0.5, that species 
with ionization potentials lower than that of hydrogen were again fully ionized, and that the 
fractional ionization of the other species was the same as that of hydrogen, i.e. 50%. 

In each case, we computed the cooling rate due to fine structure and/or metastable 
transitions from each of the listed species for a large number of temperatures in the interval 
50 < T < 10000 K and number densities in the interval 0.001 < n < 100 cm~^. In these 
calculations, we assumed that the total metallicity of the gas was 0.1 Z©, and that the gas 
was optically thin in all of the relevant transitions. We also calculated the contributions to 
the total cooling rate made by H I Lyman-a emission and Compton cooling; the rate of the 
latter was calculated for an assumed redshift z = 20. 

We used the results of our comparison to select the set of major coolants that it was 
necessary to include in our thermal model by identifying all of the coolants that contributed 
more than 25% of the total cooling rate for any of the combinations of temperature, density 
and ionization that we examined. The resulting set consisted of C, C"*", O, Si^, Compton 
cooling and Lyman-a cooling. Therefore, the only metals that we include in our chemical 
and thermal model are carbon, oxygen and silicon. The importance of carbon and oxygen 
is uns urprising - they are well known to play a major role in the cooling of the local atomic 



ISM ( IWolfire et a/.ll2003l ). and previous authors ha.ve also predicted that they will play a 
key role at high redshift (see e.g. iBronim et al.ll200ll : iBromm fc Loebll2003l : ISantoro fc Shull 



2006 



Frebel. Johnson &: Bromm 



Santoro &: Shull 



200 



3). Silicon has attracted less attention (although see 



20061 ) , but Si"*" proves to be the dominant coolant in highly ionized gas with 
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700 ^ T < 7000 K at n > 0.3 cm~^. Although gas with the very high fractional ionization 
assumed here will recombine quickly at temperatures T < 10^ K, it nevertheless seems 
prudent to include silicon in our model. Moreover, since it is necessary to include neutral 
silicon in our chemical model if we are to calculate the Si"*" abundance accurately, it requires 
little extra effort to include the effects of fine-structure cooling from Si I . 

Our choice of a 25% cutoff in this analysis is somewhat arbitrary. If we were to decrease 
the size of this cutoff, we would find that the number of species that must be included 
would increase, as both Fe"*" and S are important coolants at the 10-20% level in portions 
of our parameter space. However, as the cooling rates of the dominant coolants often have 
uncertainties that are comparable to or larger than the size of the contributions from these 
minor coolants, the accuracy we would gain by including them is less than might be expected, 
and does not (in our opinion) justify the additional complexity and computational expense 
that would be required in order to treat them. 

Finally, it is clear that our conclusions here are sensitive to our choice of elemental 
abundance ratios. For instance, an increase in the iron abundance (relative to the other 
metals) of 1-2 orders of magnitude would render Fe and Fe"*" important coolants in portions 
of our parameter space, necessitating their inclusion in our model. Such an increase would 
be expected if the enrichment of the gas were dominated by pair-in stability supernovae wit h 



masses close to the top end of the 140-260 M© allowed mass range (iHeger &: Woosleyll2002l ). 
However, studies of the abundance ratios found in extremely metal-poor stars in the Galactic 
halo, which at present give us the best picture available of the elemental composition of very 
metal-poor gas, do not fin d evidence for significant enrichrnent by pair-instability supernovae 



see e.g. the discussion in iTumlinson. Venkatesan &: Shulll l2004j ) and so at the present time 



there is no compelling reason to include iron in our chemical model. 

Our chemical model therefore consists of the eight coolants discussed above - H, H2, 
HD, C, C^, O, Si and Si"*" - together with ten additional species that play key roles in 
determining one or more of the abundances of the coolants: e~, H+, H^, H^, He, He"*", 
D, D+, O^ and Si"*"^. Our rationale for including He and He"*" in our model is that in the 
presence of a significant flux of hard UV photons. X-rays or cosmic rays, ionized helium can 
act as an important source of free electrons, and moreover can transfer charge to neutral 
carbon or silicon (but not oxygen) far more effectively than H"*" can. 

We do not include minor primordial coolants such as LiH or Hj". These are never imp or 



tant at low densities and rarely im portant at high densities (see iMizusawa. Omukai fc Nishi 



20051 . iGlover fc Savin! l2006l . 120071 ) . More significantly, we do not include molecular coolants 
such as CO or H2O, assuming instead that the bulk of the carbon and oxygen in the gas 
remains in atomic or ionized form. This assumption dramatically simplifies the chemical 
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modelling of the gas, but at the same time restricts the range of physical conditions over 
which the resulting model is useful. 

To assess the conditions for which this approximation is justified, we need to know two 
things. First, what fraction of the available carbon and oxygen must be locked up in CO 
or H2O in order for cooling from these molecules to dominate over fine structure cooling? 
Second, under what conditions are these fractions achievable within a dynamically interesting 
timescale? To answer the first of these questions, we have performed calculations using the 
treatment of carbon and oxygen fine structure cooling d iscussed in ^3.11 belqw', tog ether 



with a table-based treatment o f CO and H2O taken from iNeufeld &: Kaufman! (119931 ) and 



Neufeld. Lepp &: Melnickl ( 119951 ). We have computed /moi, defined as 



r _ Acq + Ah^o , . 

" Ac+ + Ac + Ao ' ^ ' 

where Aj is the cooling rate per unit volume due to species i, for a wide range of temperatures 
and densities. We assume that nn ^ max(nH2,ne) and that all of the cooling occurs in the 
optically thin regime. We adopt a nominal redshift z = 20 and consider only temperatures 
T > Tcmb{z) = 57.2 K, under the assumption that heating from the CMB will prevent the 
gas from cooling appreciably below this temperature. In Figure [T|, we show how /moi varies 
as a function of temperature and density. 

Figure [T] demonstrates that over most of the parameter space that we have examined, 
/mol is of order unity. We find that /moi > 10 only for temperatures very close to the CMB 
temperature (which is likely an numerical artifact, a result of the fact that our treatment of 
fine structure cooling includes the effects of radiative pumping by the CMB, while our treat- 
ment of CO and H2O cooling does not), and at high densities, where the level populations 
of the fine structure coolants start to reach their local thermodynamic equilibrium values, 
causing the fine structure cooling rate per atom to saturate. For CO or H2O cooling to be 
effective, we therefore require that about as many carbon and oxygen atoms be incorporated 
into molecules as remain in atomic form. 

With these values of /moi in hand, we can now turn to the question of whether it is 
possible to produce enough CO and H2O in the gas within an interesting timescale. Since 
we require a significant amount of oxygen to be in the form of coolant molecules, the timescale 
of interest, tconv, is given approximately by 

t ~ ^ (2) 

''conv — J-, \^ ) 

where Rf is the net rate of formation of coolants per unit volume. Therefore, to estimate 
^conv we must first estimate R{. 
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Fig. 1. — Value of /moi (see Equation [1]) as a function of temperature and density. Contours 
corresponding to /moi = 0.55 (dotted), 1.0 (dot-dashed), 3.0 (dashed) and 10.0 (solid) are 
plotted. We see that /moi < 10 for almost the whole of the parameter space considered. 



Although a complete discussion of the formation and destruction mechanisms of CO 
and H2O in low metallicity gas is beyond the scope of this paper, we will briefly summarize 
the most important points. In hot gas (T > 600 K), most CO ar id water molecules form via 



reaction pathways init iated by the hydroxyl radical, OH (see e.g. iHoUenbach fc McKeelll979 
Wagner fc Grafllll987l ). This is formed by 



O + H2 ^ OH + H, 



(3) 



but most is then destroyed by 



OH + H ^ O + H2. (4) 

However, a small fraction instead reacts to form other molecular species, such as water or 
CO, e.g. 



OH + H2 
OH + C 



H2O + H, 
CO + H. 



(5) 
(6) 



The resulting H2O and CO molecules can be destroyed by collisions with atomic hydrogen: 

(7) 



H2O 
CO 



H 
H 



OH + H2, 
C + OH. 
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If, as will generally be the case in low-metallicity gas, ^0.1 (where is the fractional 
abundance of H2 relative to the total number of hydrogen nuclei), then the destruction of 
water by reaction [7] is far more effective than its formation by reaction and water will never 
account for more than a small fraction of the available oxygen. On the other hand, the CO 
formed by reaction [6] can potentially account for almost all of the available oxygen or carbon 
- whichever is present in the smaller amount - as the destruction of CO by reaction [8] is 
ineffective at T < 5000K. Therefore, the net rate of formation of coolant molecules (primarily 
CO) is given approximately by 



^ ^mW^nom, cm-3 s-\ (9) 

where ki is the rate coefficient of reaction i, xh and xq are the fractional abundances of 
atomic hydrogen and atomic carbon, and no and n-a^ are the number densities of O and H2 
respectively. The timescale to convert significant quantities of oxygen to CO is then given 
approximately by 

t - ^ (10) 

^ 1 AgpH + ^pc ^^^^ 

For = 10~^, which is a reasonable value for low metallicity gas if H2 formation on dust 
is unimportant, this gives a timescale of approximately 

Wv^— %Myr, (12) 
where n is the number density of hydrogen n uclei, and where we have adop t ed the values for 



the various rate coefficients that are given in lLe Teuff. Millar fc Markwickl (120001 ). 



In cold gas, all of these reactions (except for reaction E]) are ineffective, and other pro- 
cesses dor ninate the formation of wa ter and CO. A good summary of the relevant chemistry 
is given in lBlack fc Dalgarnd (119771 ). When a;o+a;H2 ~ 10~^xoa;H, the most important mech- 
anism involves the formation of the OH"*" ion via the rapid ion-neutral reaction 

0+ + H2 ^ 0H+ + H. (13) 

If Xe > 1.6 X 10~^T^/^a;H2, then most of the resulting OH"*" ions simply disso datively recom- 
bine: 

OH+ + e--^0 + H. (14) 

Otherwise, they can then react further with H2 to give 1120"*" and HsO^, with dissociative 
recombination of the latter producing OH and water. CO formation follows through re- 
action [6l All of these reactions occur rapidly, and so the net rate of formation of coolant 
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molecules is given by the rate of formation of OH'*' multiplied by the fraction of OH'*' that 
is not destroyed by dissociative recombination, (1 - /dr): 

Rf - f^^^-^ - fdv)no+nu2, (15) 

For Xc = 1.6 X lO^^T^/^XHa, we have fdr = 0.5, and hence R{ = 0.5/;[j3po+nH2. In that case, 
the time required to convert most of the oxygen to OH"*" and thence to other molecules is 

^ (16) 



where we have used the fact that no+ /tiq — '^h+ / "^h owing to the rapid transfer of charge 
between oxygen and hydrogen (see §2.2.21 below). If we assume that Xe — Xii+ and that 
= 10~^, then this gives us a value for tconv of 

1000 ,^ 

icon. ^ Myr. 17 

n 

Decreasing the fractional ionization of the gas will increase tconv? but increasing it will have 
no significant effect. 

Finally, if the fractional ionization of the gas is too low for reaction [13] to operate 
effectively, then the formation of coolant molecules again occur s primarily via hydroxyl, 



which in this case is formed m ainly by direct radiative association (jjulienne. Krauss fc Donn 



197ll : ISmith fc Zweibellll976h 



O + H^OH + 7. (1^ 



In this case, tconv is simply 



^conv (.hi 



40 

~ — Gyr. (19) 
n 

Comparing the three values of tconv derived above, we see that in gas with Z = 0.1 Z©, 
tconv — lOOOn"^ Myr, regardless of the gas temperature, provided that — 10"^ and that 
the fractional ionization satisfies the constraint given above. If we compare this with the 
characteristic physical timescale of the problem of interest, tchar, then it is simple to show 
that for densities 

^ 1000 Myr ^ , , 

n < — cm-^ (20) 

t-char 

CO and H2O will not form in quantities large enough to dominate the cooling, and hence it 
is valid to ignore these molecules and all of their associated chemistry. Note also that if we 
were to include the effects of photo dissociation of OH, H2O and CO in the above analysis, 
then this would push the required density to an even larger value. 
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in a gravitationally collapsing protogalaxy, a reasonable value for tchar is the gravita- 
tional free-fall timescale tg, and in this scenario, our neglect of the molecular chemistry is 
valid as long as n < 400 cm~^. On the other hand, if we want to simulate the thermal evo- 
lution of the interstellar medium in a metal-poor dwarf galaxy, a more reasonable timescale 
may be the sound-crossing time of the disk, which is of order 100 Myr for a 1 kpc disk and 
a sound speed of 10 km s~^. In this case, our model is valid only for n < 10 cm~'^. 



As a consistency check on our conclusions here, we examined the results of lOmukai et al. 



( I2OO5I ). who model the thermal and chemical evolution of freely-falling gas at a range of 
metallicities far below solar, using a detailed treatment of the gas chemistry. They find that 
for Z = 0.01 Zq, significant conversion of carbon and oxygen to molecular form does not 
occur until n > lO^cm"^, in line with the value derived here. At lower metallicities, an even 
higher gas density is required. As another check, we have computed the evolution of the 
CO and H2O abun dances in the g as at the center of several of the simulated protogalactic 



halos discussed in IJappsen et al.l (2007; hereafter paper II), using values for the density. 



temperature, H2 abundance and abundance taken from our simulations, and modellin g 



the chemistry with the full UMIST99 chemical network (iLe Teuff. Millar fc Markwickll2000l ). 
We find that in most of these runs, our neglect of the molecular coolants is justified, as their 
abundances never become large enough for them to significantly affect the cooling. Our 
approximation begins to break down in the runs with Z = 0.1 Zq, where about 10%-20% 
of the total carbon and oxygen are incorporated into CO and H2O, which is just enough to 
affect the cooling at high temperatures and/or high densities. However, only in our runs 
with Z = Zq does it break down completely. Therefore, the use of our highly simplified 
chemical model would appear to be justified in low-density gas with Z < 0.1 Zq. 

We are thus left with a set of eighteen chemical species that must be modelled: e~, 
H+, H, H-, H+, H2, D+, D, HD, He, He+, C, C+, O, 0+, Si, Si+ and Si++. The combined 
evolution of the abundances of these species is described by a chemical network consisting 
of 74 reactions: 47 collisional gas-phase reactions (summarized in Tabled]), 12 photochem- 
ical gas-phase reactions (summarized in Table [2]), 7 grain surface reactions (summarized in 
Table [3]), and 8 reactions involving cosmic rays (summarized in Table Hj). The abundances 
of these species are also constrained by seven conservation laws: 



xh+ + + xd+ + a;He+ + a;c^ 


- + Xo+ + Xsi+ + Xsi++ 


= Xe+ xn- , 


(21) 




+ 2Xjj+ + 2xh2 + ^^HD 


= 1, 


(22) 




Xt)+ +Xy) + Xhd 




(23) 






= 3;He,tot; 


(24) 




xc + Xc+ 


= a;c,tot, 


(25) 




Xo + Xo+ 


= 2;o,tot; 


(26) 
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a^Si + a;si+ + 2a;si++ = xsi,tot (27) 

where Xi is the fractional abundance of chemical species i relative to the total abundance of 
hydrogen nuclei, and where XD.tot, ajHe.tot, a;c,tot5 a;o,tot5 and xsi^tot are the total abundances 
of deuterium, helium, carbon, oxygen and silicon, respectively. Furthermore, we assume 
in our modelling that H~ and are in chemical equilibrium, allowing us to write their 
abundances as: 

X ~ kiXuXeU (28) 

{k2X-n + k5Xii+ + ki^Xe + kiQXn + kux^+)n + R51 ' 

and 

_ {k3XKXn+ + k7Xn2X-ii+ + ki7Xu-xu+)n + CusXr^ + /^qn 
2 [kiXn + keXcju + it52 

where n is the number density of hydrogen nuclei. This assumption is generally justified 
in simulations of the cooling and gravitational collapse of gas in protogalactic halos, as the 
timescales on which H~ and reach chemical equilibrium are much sh orter than the cooli ng 



or free-fall timescale. (For a more detailed discussion of this point, see lGlover et al.ll2006l ) 



The constraints represented by equations [2TI - I29] allow us to reduce the total number of 
chemical rate equations that must be solved to only nine. In practice, we generally choose 
to solve for the abundances of the ionized species (H"*", D+, He'*', C"*", O"^, Si"^, Si"^+), H2 
and HD, but alternative choices are possible and would not significantly alter the results 
obtained. 



2.2. Selection of reactions 

The number of chemical reactions that could be included in our chemical network is very 
large, despite the limited number of chemical species involved. Fortunately, many of these 
reactions have little or no impact on the evolution of the abundances of our main coolants 
and so the number of reactions that need to be included in our chemical network remains 
reasonably small. 

We can divide the reactions that must be included into two subsets. The first subset 
consists of the reactions required to model the chemistry of hydrogen, helium and deuterium, 
including the formation and destruction of H2 and HD (reactions 1-29, 48-55, 60-63, 67-70). 
The second subset consists of the reactions required to model the carbon, oxygen and silicon 
chemistry (reactions 30-47, 56-59, 64-66, 71-74). 
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2.2.1. Hydrogen, helium and deuterium chemistry 



The amount of present in the gas is controlled by seven main reactions: collisional 
ionization of H by electrons (reaction 11), charge transfer with helium (reactions 26 & 27), 
photoionization (reaction 48), cosmic ray ionization (reaction 67), gas-phase recombination 
(reaction 13) and recombination on the surface of dust grains (reaction 61). Similar reactions 
(nos. 12, 28, 29, 49, 68, 14 & 62) partially determine the D"*" abundance, but in this case, 
charge transfer to and from hydrogen (reactions 18 & 19) is also of great importance, owing 
to the very large abundance of hydrogen relative to deuterium. Finally, the He"*" abundance 
is controlled primarily by collisional ionization (reaction 24), photoionization (reaction 50), 
cosmic ray ionization (reaction 69), gas-phase and grain-surface recombination (reactions 25 
& 63) and charge transfer with hydrogen (reactions 26 & 27). 

The remaining 24 reactions in this subset control the formation and destruction of H2 
and HD. H2 forms in the gas phase via the intermediate ions H~ and (reactions 2 & 
4), as well as on the surface of dust grains (reaction 60). It is destroyed by collisions with 
H+, e~, H and H2 (reactions 7-10), and can also be photodissociated or photoionized by UV 
radiation (reactio ns 53 & 54), or ionized by cosmic rays (reaction 70). Collisions with He 
JPove et al.lll987l ) 

H2 + He ^ H + H + He, (30) 



and He+ flBarlowlll984[ l 



H2 + He^ 



H + H+ + He, 
H+ + He, 



(31) 
(32) 



can also destroy H2, but in general these processes are not as effective as collisions with hy- 
drogen, and so they can be omitted from our simplified chemical model without significantly 
affecting its accuracy. 

The H~ and H^ ions required for gas-phase H2 formation are formed primarily by the 
radiative association of atomic hydrogen with free electrons or protons respectively (reactions 
1 & 3), and while reactions 2 & 4 generally dominate the removal of H~ and H^ from the 
gas, in hot or highly ionized gas a number of other processes become competitive (reactions 
5, 6 & 15-17) Photodissocation of H^ and H^ (reactions 51 & 52) can also become important 
if the incident radiation field is strong. 

Finally, although H P can form from intermediate ion s such as or HD+ in a manner 
analogous to H2 (see e.g. IStancil. Lepp. fc Dalgarndll998l ). most actually forms from H2 via 
reaction 20: 

H2 + D+ -> HD + H+. (33) 
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The resulting HD can be destroyed by the inverse of this reaction (reaction 21), or by 
photodissociation (reaction 55). In hot gas, the HD abundance is also influenced by reactions 
between H2 and D and HD and H (nos. 22 & 23 respectively). Note that with the exception 
of the grain surface reactions (which are discussed in §2.31 below), the reactions required to 
model the hydrogen, helium and deuterium chemistry accurately in metal-enriched gas are 
just the same as those required to model primordial gas. 

The values for many of these rate coefficients are known to within a small amount of 
uncertainty at the temperatures and densities of interest, and so our choice of the particular 
values used here should be uncontroversial. However, a few of our assumptions demand 
further comment. 

First, we note that the rates listed for reactions involving H2 or HD as a reactant 
generally assume that these molecules are not vibrationally excited. This assumption is 
reasonable at the low densities treated here, but breaks down at densities n ^ 10"^ cm~^. 

Second, we note that for several reactions involving deuterium where no readily available 
rate coefficient exists in the astrophysical literature, we have assumed that the rate is the 
same as for the analogous reaction involving H or H"*". 

The rates of processes involving cosmic rays (see Table H]) depend on the energy spectrum 
and energy density of cosmic rays. These are poorly known at the relevant energies even in 
the local ISM, and far less is known concerning high-redshift cosmic rays. For this reason, 
we do not give absolute values for these rates, but instead parameterize them in terms of 
the cosmic ray ionization rate for atomic hydrogen, (^h, which can then be considered an 
adjustable parameter of the model. 

Finally, we note that two of the reactions that regulate the H~ abundance have large 
uncertainties in their rate coefficients. The reactions in question are the associative detach- 
ment of H~ wit h H (reaction 2) ari d the mutual neutralization of H~ with H"*" (reaction 5). 



As discussed in lGlover et al.l (120061 ) . the uncertainties in the rates of both of these reactions 
may be as large as an order of magnitude. In gas with a high fractional ionization, this 
uncertainty can lead to a significant uncertainty in the H2 formation rate and in the final 
H2 fractional abundance, particularly in the presence of a strong ultra violet background ra - 



diation field. In the chemical model presented here, we ha ve followed iGal 



and have adopted a rate coefficient f or reaction 2 taken fro m iLaunay et al.l ( 



coefficient for reaction 5 taken from iMoseley et al.l (jl970l ). However, we caution the reader 



i fc Fallal fll998h 
I991I ) and a rate 



that this should not be regarded as an endorsement of the accuracy of these particular values. 
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2.2.2. Carbon, oxygen and silicon chemistry 

As discussed in §2. II above, we do not include the molecular chemistry of these elements 
in this simplified model. Our treatment of the carbon, oxygen and silicon chemistry is 
therefore purely a treatment of the charge balance of these species. 

We begin with oxygen, in many respects the simplest of the three to treat. The ionization 
potential of neutral oxygen is only 0.02 eV larger than that of neutral hydrogen, and so 
charge transfer between O"*" and H or H"*" and O (reactions 36 & 37) occurs rapidly. Since 
the hydrogen abundance is orders of magnitude larger than the oxygen abundance, this 
means that the ratio of ionized to neutral oxygen is controlled by the ratio of ionized to 
neutral hydrogen, i.e. that 

I • 34 

In most circumstances, reactions 36 & 37 are the only reactions required in order to accurately 
model the oxygen chemistry. Nevertheless, for completeness we also include several other 
processes: radiative recombination (reaction 32), grain surface recombination (reaction 65), 
coUisional ionization (reaction 35), charge transfer with He^ (reaction 38), photoionization 
(reaction 57) and cosmic ray ionization (reaction 72). 

In the case of carbon, the situation is rather different. First, charge transfer from H'*' 
to C (reaction 39) is much less effective than charge transfer from H"*" to O, and so reaction 
39 plays a far less important role in the carbon chemistry than reaction 37 does in the 
oxygen chemistry. Second, carbon has an ionization potential of only 11.26 eV. This means 
that charge transfer from C+ to H is significantly endothermic, rendering it unimportant 
at T < 10^ K (although we include it here for completeness). It also means that neutral 
carbon can be photoionized by ultraviolet photons with wavelengths A > 912 A, which can 
penetrate easily into low metallicity protogalactic gas. Consequently, photoionization of C 
(reaction 56) plays an important role in the carbon chemistry, whereas photoionization of O"*" 
is unimportant outside of H ll regions. In addition to charge transfer and photoionization, 
ionized carbon can also be produced by coUisional ionization (reaction 33), although this 
is important only for temperatures T > 9000 K, by charge transfer with He'*' (reaction 
41) and by cosmic ray ionization (reaction 71). C"*" is removed from the gas primarily by 
recombination in the gas phase (reaction 30) and on grain surfaces (reaction 64), although 
if the free electron abundance is small, charge transfer from C"*" to Si (reaction 44) can also 
become important. 

In the case of silicon, a similar set of processes operate: coUisional ionization (reaction 
34), cosmic ray ionization (reaction 73), gas-phase recombination (reaction 31), grain surface 
recombination (reaction 66), charge transfer from H"^, He^ or C+ to Si (reactions 42, 43 & 
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44), and photoionization by an ultraviolet radiation field (reaction 58), assuming one is 
present. In the case of silicon, however, we also include the doubly ionized ion. Si"'""'", in our 
chemical model, since this can be produced by charge transfer from H"*" to Si"'" (reaction 45) 
with an endothermicity of only a few eV, rendering it potentially important at temperatures 
~ 10^ K. It can also be produced by photoionization (reaction 59), although in this case the 
16.35 eV energy requirement renders this process important only within H II regions, or by 
cosmic ray ionization (reaction 74). The Si^"'" produced by these processes can be destroyed 
by charge transfer with neutral hydrogen (reaction 46), or by recombination (reaction 47). 

Regarding the accuracy of our adopted reaction rates, we note that while many have 
well determined rate coefficients, there are two notable exceptions. The only available rate 
coefficients for reactions 43 and 44, which involve the transfer of charge from He"'" to Si and 
from C+ to Si respectively, assume that the reactions both proceed at the Langevin rate. In 
practice, highly exothermic charge transfer reactions often proceed at a rate far below the 
simple Langevin rate, and so the rates of these reactions may be overestimated in our model. 



2.3. Grain surface chemistry 



If dust grains are present in the metal-enriched gas, then our chemical model must 
account for the effects of reactions occurring on the surface of the grains, as these reactions 
are often far more effective than their gas-phase counterparts. We therefore include a small 
number of grain surface reactions in our model, as summarized in Table [31 

The only grain surface reaction between neutral species that is included in our model 
is the formation of molecular hydrogen (reaction 60). This is a hugely important reaction 
in the local interstellar medium, since at typical interstellar densities it is the only process 
capable of producing t he large quantities of molecular hydrogen that are directly observed 
(see e.g. 



Wakkeil l2006h or inferred from other molecular tracers, such as CO. Theoretical 



modeling (IHirashita fc Ferraral l2002l : iGloverl l2003l : ICazaux &: Spaand |2004| ) suggests that it 
remains important down to metallicities of order 10~^ Z©, within the range of applicability 
for the chemical model presented here. 



The rate that we adopt for this process is based on the widely used rate of lHollenbach fc McKee 



( I1979I ). This was derived for solar m e tallici ty gas, assuming a distribution of grain sizes as 



given in iMathis. Rumpl fc Nordsieckl (119771 ). and to adapt it for use in low metallicity gas. 



we simply assume that the rate scales linearly with the metallicity Z. The validity of this 
assumption is open to question, as many of the features of the grain population, such as the 
grain size distribution or the mix of compositions, may differ greatly between Milky Way dust 
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and protogalactic dust. An alternative, physically motivated approach would be to adopt a 
grain size distribution and mix of compositions based on the results of numerical modeling 



of dust formation in high redshift su pernovae (iTodini fc Ferraral l200ll : iNozawa et al.l 12003 



Schneider. Ferrara &: Salvaterrall2004j ) and then to cor apute the H2 forniation rate expected 
for this grain population. This is the approach used by lSchneider et al.l (120061 ). However, the 
uncertainties associated with this approach are considerable. To begin with, the predictions 
of the numerical mod els are highly sensiti ve to the degree of mixing assumed to occur within 
the supernova ejecta (INozawa et al.ll2003l ). In addition, dust destruction in the reverse shock 
is typically not taken into account in these models, and to the best of our knowledge, the 
amount of dust that survives the passage of the shock has yet to be fully quantified (although 
research in this area is actively proceeding; R. Schneider, priv. comm.). Finally, the changes 
wrought on the grain popula tion by subsequent processing i n the interstellar medium or 
intergalactic medium (see e.g. IVenkatesan. Nath &: Shullll2006l ) are not fully understood. In 
light of these uncertainties, we do not believe that this approach currently offers much of an 
advantage over our simple assumption of an H2 formation rate that scales with metallicity. 

We do not include any neutral-neutral surface reactions in our model other than H2 
formation. For neutral-neutral surface reactions to be able to significantly affect the ability 
of the gas to cool, they must be able to alter the abundances of neutral carbon, oxygen or 
silicon by a substantial amount. We can place an upper limit on the rate at which these 
reactions occur if we assume that the reaction probability /react = 1, i-e. that every collision 
between a metal atom and a grain results in a reaction. In that case, the reaction rate per 
unit volume for an atomic species i is given by 



-Ri = VthiArii, 



(35) 



where Vth,i is the thermal velocity of atomic species i, and A is the tota l surface area of grains 
per un it volume of gas. For Milky Way dust of the type assumed by iHoUenbach fc McKee 
( 119791 ). y4 ~ 3 X 10~^^r;, cm~^, where n is the number density of hydrogen nuclei. At lower 
metallicity, our assumption that the grain size distribution does not change with metallicity 
implies that the value of A in gas with a metallicity Z is simply A ~ 3 x lO~^^(Z/Z0)ncm~^. 
We can therefore rewrite equation [321 as 



Ri ~ 4.7 X 10 
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T 
rrii 



1/2 



Z 



nrii, 



(36) 



where is the mass of species i in atomic mass units, 
timescale tconv = ni/Ri is then 



The corresponding conversion 



2.1 X 10^%-^ ( — 



-1/2 



Z 



-1 



(37) 
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If /react < 1, then this expression becomes: 



Wv = 2.1 X 10^V-a>"' ' s. (38) 



^0 / 

If we equate this to a characteristic physical timescale tchar, then we can show that tconv > ^char 
as long as 

For gas with a temperature T = 2000 K and metallicity Z = 0.1 Zq, and with atomic carbon 
as the colliding species, this gives: 

500/-^^,tMyr 3 
n < — — cm'^. (40) 

^char 

At higher temperatures, the limiting n will be slightly smaller, but at T > 2000 K, collisional 
dissociation of most molecular species is highly effective, and so in this temperature regime, 
grain surface reactions are unlikely to be important. 

From this analysis, we see that even if /react = 1, neutral- neutral grain surface reactions 
are unimportant at gas densities n < (500 Myr/tchar)(Z/0.1 Zq)~^ cm~^. In gravitationally 
collapsing gas with t^har = this corresponds to n< lOO(Z/O.lZ0)~^cm~^. Note that a key 
point here is that for neutral-neutral grain surface reactions involving our atomic coolants 
(C, O etc.) to significantly affect the thermal behaviour of the gas, a large fraction of the 
total number of coolant atoms must react, whereas for grain surface reactions to affect the 
H2 cooling rate by significantly altering the H2 abundance, only a small fraction of the total 
number of hydrogen atoms must react. This means that in the case of H2, the relevant 
conversion timescale is several orders of magnitude shorter, and hence our density limit 
would be significantly smaller if we were to omit grain surface H2 formation from our model. 

One way in which this simple analysis could break down is if /react were much larger for 
reactions involving the coolant atoms than for H2 formation. However, a recent analysis of H2 



formation on grain surfaces by lCazaux fc Tielend (120041 ) that takes both physisorbed (i.e. van 



der Waals bonded) and chemisorbed (i.e. chemically bonded) hydrogen into account demon- 
strates that in the conditions of interest in this paper, /react ~ 1 (wi th the result that the 



compu ted H2 formation rate is very similar to the widely-used rate of iHoUenbach fc McKee 



(I1979I ) that is used in our model). A comparable analysis has not been performed for reac- 
tions involving C, O or Si, but clearly they cannot have /react > 1, and so there is little scope 
for these reactions to occur significantly faster than H2 formation. 

Finally, we include in our surface chemistry model six important reactions involving 
ions: the recombination of II+, D+, He^, C^, and Si"^ with electrons on the surface of 
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grains. As grain surface recombination is a non-radiative process, it proceeds at a much 
faster rate than radiative recombination in the gas phase. Moreover, if a typical grain is 
negatively charged, then the effective cross-section for collisions will be much enhanced over 
the geometric cros s-section due to Coloum b focussing. Grain charging is largely determined 
by the parameter (iBakes fc Tielenslll994l ) 



, G^JT 
%l) = , 



rip 



(41) 



where G ~ O. OlJ^i is a mea sure of the radiation energy density between 6 eV and 13.6 eV 
relative to the lHabing] (119681 ) field. When ip is small, most grains are negatively charged, and 
so in these conditions grain surface recombination can be important even if neutral-neutral 
grain surface reactions are unimportant. 



To model grain surface recombination, we follow IWeingartner fc Draind (l2001al ). We 
adopt their rate coefficients for the recombination of H"*", He"*", C"*" and Si"*". For O"*" and D"^, 
we use the facts that O"^ and D"*" have almost the same ionization potential as H"*' and that 

— 1/2 

the ion arrival rate at the grain scales as m- (where rrii is the ma ss of the ion in atomic 
mass units) to derive rates by appropriately scaling the H"*" rate. The IWeingartner fc Draine 
rates were all computed for Milky Way dust, and to adapt them for use in low metallicity 
gas, we again assume that they scale linearly with Z, with the same caveats as before. 



2.4. Photochemical rates 

In Table |2l we list the cross-sections for all but two of the photochemical reactions 
included in our model. The two exceptions are H2 photodissociation and HD photodissocia- 
tion, which are caused by absorption in a large number of discrete spectral lines, and which 
are discussed separately in §2.4. II below. 

Given the cross-section, (y{E), the corresponding photochemical rate can be obtained 

from 

i?photo = 47r r ^i^)lMe-riE) [1 + f^E)] dE (42) 

Jeo ^ 

where E is the photon energy, Eq is the energy threshold, I{E) and t{E) are the mean 
intensity (in units of eV s~^ cm~^ eV~^ sr~^) and optical depth for a photon of energy E 
(both of which are problem dependent), and where f{E) is a factor that accounts for the 
effects of secondary ionizations. It is generally a reasonable approximation to se t f(E) = 



unles s the ionizing compon ent of the radiation field is dominated by X-rays (lAbel et al. 



19971 : iGlover fc Brand! l2003l ). In the case that X-rays dominate, fits for f{E) for H and 
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He photoionization as a fun ction of the fractional ionization of the gas can be found in 
ShuU fc van Steenberg (Il985l ) and iDalgarno. Yan &: Liul (Il999l ). The effects of secondary 
ionizations on the other processes hsted here are generally negligible, owing to the large 
abundances of neutral H and He relative to all other species. 

If the gas is optically thick at the He photoionization threshold, then an additional 
process that must be taken into account is the photoionization of H by the diffuse emis- 
sion produced by He"*" recombination. In the limit of high optical depth, the on-the-spot 
approximation ap plies, and we can model this process as a local ionization rate with a value 
JOsterbrocklll989h 



-Rpi = [2//^25,rr,A + (0.96 - ?/)/C25,rr,B + fc25,di] nen^^+ Cm ^ S' 



where y is given by 



y 



(43) 



(44) 



'^HO"48(-E'th,He) + '"'HeO"5o(-E'th,He) ' 

where -Eth,He = 24.6 eV is the He ionization threshold. (Note that if n-ne/nn ~ 0.08, as is the 
case in primordial gas with a low fractional ionization and low H2 abundance, then y ~ 0.68). 
If the on-the-spot approximation does not apply, then the radiative transfer of this diffuse 
emission must be modelled in some fashion. However, a discussion of appropriate techniques 
for doing so lies well beyond the scope of this paper. 

In some circumstances it may also be necessary to take account of the photodissociation 
of H~ and H^ by photons produced by ionized hydrogen and helium (both recombination 
emission and bremsstrahlung) . This process is generally important only in neutral gas close 
to a significant volume of dense, ionized gas (e.g. in the neu tral gas immediately surrounding 
an expanding ionization front). It is discussed in detail in iGloverl (120071 ) and so we do not 
discuss it further here. 

Finally, we note that for each of the photoionization or photodissociation rates listed in 
Table [2l there is a corresponding photoheating rate, given by 



R 



■heat 



Eo 



a{E)I{E) 
E 



e— (^)(^ _ Eo)7]{E - Eo)dE. 



(45) 



where E — Eq is the energy of the primary photoelectron and ri[E — Eq) < 1 gives the 
fraction of this energy conve rte d to heat, which can be ca lculated using the results of 
ShuU fc van Steenberg) (119851 ) or IDalgarno. Yan fc Liul (119991 ). In practice, photoheating 
from the photoionization of H and He usually dominates over the other contributions by a 
wide margin. 
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2.4- 1- H2 and HD photodissociation 



Although the binding energy of H2 is only 4.48 eV, photons of this energy are not able to 
dissociate H2 effectively, as the simplest dissociative tran sition - excitation to the vibrationa l 
continuum of the ground state - is strongly forbidden (IField. Somerville fc Dresslerlll966l ). 
Transitions to the repulsive b'^S^ state, the least energetic of the excited electronic states of 
H2, are also forbidden, and so photodissociation takes place primarily through excitation to 
the Lyman (B^S+) or Werner (C^n„) electronic states followed by radiative decay to the 
vibrational continuum of the ground state. As a number of vibrational levels are accessible 
in each excited state, photodissociation takes place thr ough a number of discrete absorption 
lines, known as the Lyman and Werner band systems (jStecher fc Williams! 119671 ). 



In optically thin gas, the photodissociation rate can be written as 



diss 



■diss,v,J fv,J 



(46) 



where f^^j is the fraction of H2 molecules that have vibrational and rotational quantum 
numbers {v, J) in the electronic ground state, and Rdiss,v,j is the photodissociation rate due 
to transitions out of {v, J). The latter can be written as 



(47) 



v',J' 



where Cv,j,v',j' is the pumping rate from level {v, J) in the electronic ground state to level 
(f', J') in either the Lyman or Werner states, fdiss,v',j' is the fraction of decays from v', J' 
which end in the vibrational continuum of the ground state (rather than back in some 
bound state), and where we sum over all accessible levels. Given appropriate molecular 
data, calculation of Rdiss,v,j is straightforward for each bound level {v, J) in the electronic 
ground state. To then calculate -Rdiss, one also needs to know the level populations f^j. 

If we assume that the mean intensity /(z/) = hI{E) is independent of energy, and that 
all of the H2 is in the f = 0, J = level (i.e. the para-hydrogen ground state), then _Rdiss 
evaluates to 

i?diss = 1-38 X 10^ I {u) s-\ (48) 

where we have made use of molecular data taken from I Ab grail et al. Jl993al lb[) and lAbgrall. Roueff fc Driral 
(I2OOOI ). This expression remains a good approximation in the more general case that /(i^) is 
allowed to vary with frequency, provided that the variations ar e not too extreni e and that 
we replace /(i/) in equation HSl with /(i^), where = 12.87 eV (lAbel et al.l 119971 ): 



i?diss = 1.38 X lO^I(u) s-\ 



(49) 
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RelaxiiiiB; the assumption that all of the H2 has J = also makes little differ ence to -Rdi ss 
( lGloverll200ll ). Vibrational excitation of the H2 makes a much larger difference (IShulllll978l ). 
but at low gas densities we would expect the populations of the vibrational levels of H2 to 
be very small. Equation SHI therefore gives a reasonable estimate of the optically thin H2 
photodissociation rate within the regions of parameter space for which our chemical model 
is valid. 

If enough H2 is present in the gas, then the Lyman- Werner lines can become optically 
thick, leading to a reduction in the H2 photodissociation rate, an effect known as H2 self- 
shielding. If the gas is at rest, th en the effects of H2 self-s hielding can be treated quite 
accurately using the prescription of iDraine &: Bertoldil (119961 ). They parameterize the self- 
shielding with a shielding function fsh, defined to be the ratio of the H2 photodissociation 
rate in self-shielded gas to the rate in optically thin gas. They demonstrate how to calculate 
fsh as a function of the gas temperature and the H2 column density and also construct the 
following useful fitting function: 



0.965 



sh 



+ 



0.035 



■exp [-8.5 X 10-^(1 + 



(l + s)i/2 
^, A^H2 is the H2 column density, 



(50) 



where x = Nii^/b x 10 cm~ , A^H2 is the H2 column density, = 6/10 cm s~ and b 
is the Doppler broadening parameter. Although iDraine &: Bertoldil (119961 ) assume a semi- 
inifinite slab geometry in their models, their approach is easy to extend to more complicated 
geometries. 

Unfortunately, their simple treatment breaks down in gas which is not at rest. Doppler 
shifts due to the motions of the gas cause H2 in different regions to absorb at slightly 
different wavelengths, and if these Doppler shifts are comparable to or larger than the thermal 
linewidth of the gas (as will be the case in transonic or supersonic gas respectively), then the 
effect is to reduce the amount of self-shielding that occurs. An accurate treatment of H2 self- 
shielding in this regime probably requires one to solve the full frequency-dependent transport 
equation, which cannot currently be done in a computationally efficient manner within a 
three-dimensional hydrodynamics code. Consequently, various different approximations have 
been used to study H2 photodissociation in this regime. 



The simplest approach is to ignore self- shielding entirely (see e.g. iMachacek et al.ll2001 



2OO3I ). This is a good approximation if the velocities in the gas are large and the H2 column 
densities are small, but otherwise will significantly overestimate the ph otodissociation rate . 



At the other extreme, one c an ignore the effects of Doppler shifts (e.g. lYoshida et al.l 12003 



Hosokawa fc Inutsukal l2006l ). This is a good approximation if the H2 column density is 
sufficiently large (A^H2 > lO^^cm"^) that the Lorentz wings of the Lyman- Werner line profiles 
dominate the line widths, as in this case the line widths will be much larger than any likely 
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Doppler shifts within the molecular gas. On the other hand, this approach will underestimate 
the true photodissociatii 
low-metallicity systems. 



the true photodissociation rate when A^H2 < 10^^ cm ^, as is the case in many interesting 



Another approximation has recently been suggested by lAhn fc Shapiro! (120071 ) . They 
use equation (1501) to compute /sh, but adopt a value for h that includes both a thermal 
contribution and one arising due to the velocity dispersion of the gas. In practice, this 
means that they treat their H2 as having an effective h equivalent to that in a purely thermal 
gas with T = lO'^ K. The accuracy of this approximation depends on the correlation length 
of the velocity field. If this is small compared to the other length scales of interest, then 
treating the velocity dispersion in this fashion is reasonable and should give a fairly accurate 
result. On the other hand, if the velocity field is dominated by large-scale bulk motions (such 
as infall into a protogalaxy) , then this approximation will be significantly less accurate. 



Finally, in iGlover et al.l (120061 ) and in paper II, we use a local approximation in which 
only the H2 within a single SPH srn oothing length is assumed to contribute to the shielding 
see also Iciover k Mac LowlboOTal lbl. for a grid-based version of this approach). This fairly 
crude approximation is intended to take account of the fact that H2 close to a given point of 
interest is more likely to have only a small relative velocity than gas a large distance away. 
It will generally underestimate the amount of self-shielding, but nevertheless represents an 
improvement over neglecting self-shielding entirely. Aside from its inevitable inaccuracy, 
this approximation also suffers from the disadvantage of being resolution dependent, as 
increasing the number of SPH particles in the simulation will generally decrease all of the 
SPH smoothing lengths and hence will cause a systematic increase in f^h- On the other 
hand, it has the significant advantages of being computationally efficient (as only local data 
is required) as well as being very easy to implement. 

To sum up, a number of different approximate methods exist for treating H2 self- 
shielding in large numerical simulations, but none are entirely satisfying. Further work 
on this problem is definitely called for. 

Turning now to HD, we note that HD photodissociation in optically thin gas can be 
treated in much the same way as H2 photo dissociation. The necessary molecular data for 
HD can be found in lAbgrall fc Rouefll (120061 ). and the resulting photodissociation rate for a 
radiation field with a flat spectrum can be written as 



diss,HD 



1.5 X Kfliv) S 



(51) 



10i3cni-2, 



which is only ~ 10% larger than the H2 rate. Above an HD column density A^hd 
self-shielding of the HD lines significantly reduces the ph otodissociation rate . For a static gas, 
this process can again be modeled using the approach of iDraine fc Bertoldil (Il996l ). although 
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the same problems arise when one tries to extend this approach to a gas distribution which 
is not static. However, in the case of HD, we face an additional complication: if the HD 
column density is sufficiently high for HD self-shielding to be significant, then the H2 column 
density will be very much larger (since even in significantly fractionated regions, one typically 
has an HD:H2 ratio of no more than about 10"'^). Consequently, the line widths of the H2 
Lyman- Werner lines are not negligible, and some degree of overlap between these lines and 
the HD absorption lines will occur. Additionally, if the H2:H ratio is small, as will often 
be the case in the systems of interest, then a significant HD column density implies a large 
neutral hydrogen column density, which means that absorption of radiation in the Lyman 
series lines of atomic hydrogen must also be taken into account. 

These effects are difficult to include accurately in a simple treatment of HD self-shielding 
and before attempting to do so it is reasonable to ask whether an accurate treatment of 
HD self-shielding is really required. We argue that in many cases of interest it is not. 
Comparison of the rate at which HD is photodissociated in optically thin gas with the rate 
at which it is destroyed by reaction 21 demonstrates that the latter dominates whenever 
nH+ ^ 10~^J2i('^)7 where J2i(z^) is the strength of the radiation field at hi? = 12.87 eV in 
units of 10~^^ ergs s~^ cm~^ Hz~^ sr~^. Photodissociation therefore dominates only when the 
UV field is strong or the proton number density is small. However, in either case, it is difficult 
to see how the large column densities of HD and H2 required for effective shielding could be 
built up or maintained. Therefore, we suspect that for most applications, treating the HD in 
the optically thin limit is probably sufficient, as in the conditions where this approximation 
breaks down, photodissociation is unlikely to be important. 



3. Thermal processes 

3.1. Fine structure cooling 

As we do not include molecular coolants such as CO or H2O in our chemical model of 
metal-enriched gas, for the reasons outlined in §2.11 the main contribution that the metals 
make to the cooling of the gas is through fine structure line emission from neutral C, O and Si 
atoms and C"*" and Si"*" ions. To model this emission, we assume that the populations of all the 
electronically excited levels of these atoms and ions are negligible, an approximation which 
should be highly accurate at the gas densities considered in this study. This assumption 
allows us to model C"*" and Si+ as two-level systems and C, O and Si as three-level systems, 
allowing us to compute their effects in a straightforward fashion. For a two-level ion, if we 
denote the ground state as level and the excited state as level 1, then the power radiated 
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per unit volume can be written as 

A = (Aio + 5io/io)^iorii, (52) 

where ni is the number density of ions in level 1, Aiq is the Einstein coefficient for spontaneous 
emission for the transition from level 1 to level 0, Biq is the corresponding coefficient for 
stimulated emission, Eiq is the energy of the transition, and Jio is the mean specific intensity 
at the frequency of the transition. If Iiq ^ 0, then the ions will also absorb energy from the 
radiation field, at a rate 

r = BoJioEiono, (53) 

where uq is the number density of ions in level and Bqi is the Einstein coefficient for 
absorption from level to level 1, which is related to Biq by Bqi = ( (71/(70) -Bio, where go and 
gi are the statistical weights of levels and 1 respectively. The net loss of energy per unit 
time per unit volume is therefore 

A' = A-r 

= Eio ({^10^1 + BiqIiqItii - igi/go)no]}) , (54) 

and it will be seen that if 

Aio + BiqIiq 
{gi/go)BioIio 

then the ions will absorb more energy than they emit, and so the gas will actually gain 
energy. 

To compute A', we need to know several pieces of atomic data - the values of Aiq, Eiq, 
go and gi, which are summarized for C"'' and Si"^ in Table [5]- together with the values of Iqi, 
hq and rii. To compute uq and ni, we assume that the levels are in statistical equilibrium, 
in which case: 

(fioi/01 + CQi)nQ = {Aio + BiqIiq + Cio)ni, (56) 

where Cqi and Cio are the total rates of collisional excitation and de-excitation respectively. 
These are related by 

and so once one is known, the other can be computed easily. In Table [6], we list collisional 
de-excitation rates for collisions between C"*" or Si^ and various possible collision partners 
such as H, H2 or e~. Given the number densities of these species, Cio can be easily computed, 
since 

Cio = ^gicfc^fc, (58) 
it 
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where qio,k is the colhsional de-excitation rate for a colhsion with chemical species k with 
number density n^. For C^, we include the effects of collisions with electrons, atomic hy- 
drogen and molecular hydrogen (in both ortho and para forms). The collision rate with H"*" 
is negligible at the temperatures of interest due to the strong Coloumb repulsion, and the 
abundances of the other chemical species included in our model are too small for them to be 
important collision partners. For Si"*", we include only the effects of collisions with H and e~, 
as rates for collisions with H2 are not available. However, provided that the H2 abundance 
is small compared to the atomic hydrogen abundance, this is unlikely to be a major source 
of error. 

Finally, to compute Jio, we assume that the only significant radiation field present at 
the infra-red and sub-millimeter wavelengths of the fine structure transitions is the cosmic 
microwave background. In that case, Jio is simply given by the value of the Planck function 
at the frequency of the transition for a radiation field with temperature T = Tqmb = 
2.726(1 + 2) K. 

For the three-level atoms (C, O and Si), we use a very similar approach. In this case, 
the power radiated per unit volume is 

A = (Aio + 5io/io)^io^i + (^20 + ^20/20)^20^2 + (^21 + B2ihi)E2in2, (59) 

and the power absorbed per unit volume is 

F = BqiIiqEiqUq + 502/20 -E'20'^o + Bi2l2iE2ini, (60) 

where 0, 1 and 2 denote the ground state and the two excited states respectively. The level 
populations no, rii and n2 are found by solving 

(Soi/01 + Cm + S02/02 + Co2)no = (Aio + 5io/io + Cio)ni + (A20 + ^20/20 + C^2o)n:(pl) 
(5io/io + Cio + 5i2/i2 + Ci2)ni = (5oi/oi + Coi)no + {A21 + 521/21 + ^^21)^2, (62) 

with all symbols having their obvious meanings. To compute the total collisional excitation 
and de-excitation rates for carbon and oxygen we include the effects of collisions with ortho 
and para-H2, atomic hydrogen, protons and electrons. For silicon, we include only the effects 
of collisions with H and H"*", as rates for collisions with H2 or electrons do not appear to be 
available. The rates used are summarized in Table El 

3.2. Other coolants 

Apart from fine structure emission, we also include in our thermal model several other 
processes that can lead to the cooling of the gas. These are summarized in Table [3, along 
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with a reference to the source (or sources) from which the associated coohng rate has been 
taken. In most cases, the rate itself is also listed. 

In hot, ionized gas, cooling is dominated by electron impact excitation of atomic hydro- 
gen (Lyman-a cooling), atomic helium and He"*". Excitation of atomic helium occurs from 
both the l^S ground state and the 2'^S metastable state. In common with previous authors, 
we assume that the population of the 2^S state is set by the balance between radiative recom- 
bination to triplet states and radiative decay to the ground state. One consequence of this 
assumption is that the number density of He atoms in the 2^S state, 'n.He{23s)! is proportional 
to the product of the number densities of free electrons and of He'*', i.e. 'n.He(23s) oc nenHe+, 
which means that the cooling rate from metastable helium scales as nHe(23s)'^e_2L,'^Ilhie+- 
To model cooling from H, the He metastable stat e, and He"*", w e use rates from ICenI (119921 ). 
which themselves were based on earlier work by Blackl (1981 ). To m odel cooling from the 
He ground state, we use our own fit to the data of iBray et al.l (120001 ). 



A number of other processes are of importance in ionized gas. We include cooling due 
to collisional ionization of atomic hydrogen and atomic helium, the gas phase recombination 
of H"*" and He"*", ionic recombination on dust grains, Compton scattering of CMB photons by 
free electrons (Compton cooling), and thermal bremsstrahlung. The rates adopted for all of 
these processes are summarized in Table [71 As the chemical and thermal model presented 
here is not designed to be used for the study of very hot gas, we do not inc lude processes 
involving He"'"^ (although rates for these processes can be found in lCerull992l ). 



In neutral gas, all of the aforementioned processes become ineffective. In cool, neutral 
gas, most of the cooling comes from H2 or from the fine structure lines of carbon, oxygen and 
silicon, which have already been discu ssed above. Cooling fro m H2 is treated in our model 
through use of the cooling function of iLe Bourlot et al.l (119991 ) . which we have extended to 
temperatures below 100 K by assuming tha t only the J = 2 — > and J = 3 — > 1 transitions 



contribute significantly to the cooling rate. iLe Bourlot et al.l (119991 ) tabulate the H2 cooling 
rate as a function of temperature, density, H:H2 ratio and ortho:para ratio. For simplicity, in 
our implementation we do not track the evolution of the H2 ortho:para ratio, instead keeping 
it fixed at 3:1, but we note that variations in this ratio are unlikely to significantly affect the 
H2 cooling rate at temperatu res at which it contribu tes significantly to the total cooling rate 



see, for instance, figure 5 in lLe Bourlot et al.lll999l ). A comparison of the lLe Bourlot et al. 



H2 cooling rate with various other rates that have been used in the literature is given in 
paper II (Figure 1). 

At low temperatures (T< 200 K), HD cooling becomes increa singly important an d can 
dominate the total cooling rate if sufficient fractiona tion occurs (IGalli fc Fallal 119981). To 
model HD cooling, we use the recent cooling function of lLipovka. Nunez-Lopez, fc Avila-Reese 
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( 120051 ) ■ They provide a complicated fit as a value of temperature and density that is valid 
over a wide range of both. For further details, the interested reader should consult their 
paper. 

We include two final processes that can become important in some circumstances. Cool- 
ing due to H2 collisional dissociation is modelled under the assumption that each dissociation 
removes ~4.48 eV of thermal energy from the gas. It is an effective source of cooling only at 
temperatures above a few thousand Kelvin. In practice, the H2 abundance in zero metallicity 
or low metallicity gas is frequently too small for this process to be important. 

We also incl ude the effects of energy trans fer from the gas to the dust grains (if present), 
using a rate from iHoUenbach &: McKed (119891 ). At solar metallicity, this becomes important 
at a density n ~ lO^cm"^, but at lo wer metallicities, it d oes not dominate until significantly 
larger densities are reached (see e.g. lOmukai et al.l 120051 ). 



3.3. Heating 



We include in our model the effects of several processes that can heat the gas. Most 
of these processes operate only if a radiation background is present. The first of these is 
photoelectric emission from dust grains. This operates as follows: photons that interact with 
dust grains can cause the ejection of energetic electrons from the grain if the photon energy 
exceeds the work function of the grain. As the energy carried by the electrons is quickly 
thermalized, this leads to the heating of the gas. This process is of great ir nportance in the 
local ISM and has been examined in detail by a number of authors (e .g. iBakes fc Tielens 
19941 : IWolfire et allll995l : IWeingartner fc Draineir2001bl : iRae et al.l[20ol . 



To accurately compute the effects of photoelectric emission, we need to know the grain 
size distribution and the composition of the grains. However, as discussed previously, large 
uncertainties exist concerning the properties of grains in low metallicity protogalactic gas. 
We therefore make the same assumption here as we did in our treatment of grain surface 
chemistry, i.e. that the dust has the same properties as Milky Way dust, but has an abun- 
dance that is reduced by a factor (Z/Z0). This assump tion allows us to us e the following 
expression for the photoelectric heating rate, taken from lWolfire et al\ (119951 ): 



pe 



1.3 X 10 



-24, 



G 



Z 



n erg s ^ cm ^ , 



(63) 



where G ~ 0. 01J:;.i is a mea sure of the radiation energy density between 6 eV and 13.6 eV 
relative to the iHabina (119681 ) field, and where e is the photoelectric heating efficiency, given 
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by 



4.9 X 10" 



+ 



3.7 X lO-^TO-^ 



1.0 + 4.0x10-3^0-73 1.0 + 2.0 X 10-V' 



(64) 



where T4 = T/10^ K. The parameter ijj that c ontrols the photoelectric heating efficiency is 
given by ip = G^/T/0.5ne (IWolfire et a/.ll2003l ): note that this differs by a factor of 2 from 
the parameter ip introduced in §2.3[ For small ip, most grains are negatively charged and 
e ~ 4.9 X 10"^ + 3.7 x lO-^T^-^. On the other hand, for large ip, most grains are positively 
charged, and e is small, as it is difficult for photons with energies E < 13.6 eV to dislodge 
further electrons from the grains. 

A second source of radiative h eating is the photodissoc iation of H2. Our treatment 
of photodissociation heating follows iBlack fc Dalgarnol (119771 ): we assume that each pho- 
todissociation deposits 0.4 eV of themal energy into the gas. As well as photodissociat- 
ing some of the H2, an ultraviolet background will also produce vibrationally excited H2 
via radiative pumping of the excited levels. In dense gas, this pumping leads to heating 
as most of the excited molecules undergo coUisional de-excitation. We include the effects 
of radiative pumpin g by adopting a pumpin g rate that is 8.5 times larger than the pho- 
todissociation rate (jPraine &: Bertoldil Il996l ). and assuming that each excitation transfers 
an average of 2 (1 + ncr/n)~^ eV to the gas (iBurton et al.l Il990l ). where is the crit- 
ical density at which collisional de-excitation of vibrationally excited H2 occurs at the 
same rate as radiative de-excitation. Our value for n„ is a weighted harmonic mean of 
the value f or H2-H collisions given by iLepp &: Shulll (Il983l ). reduced by a factor of ten as 
advised by Martin. Schwarz fc Mandy ( 19961 ). and the value for H2-H2 collisions given by 
Shapiro fc Kanej Jl987h . 



Heating due to the photoionization of H or He has already been discussed in §2.41 and 
we do not discuss it further here. As is also discussed in that section, we do not include 
heating due to the photodetachment of H", photodissociation of H^ or photoionization of 
C, O, Si or Si"*", as the contribution from these processes is not significant. 

We also include the effects of heating due to H2 formation. The formation of an H2 
molecule via reaction 2 releases 3.53 eV of energy, while formation via reaction 4 releases 
1.83 eV, and formation on a grain surface (reaction 60) releases 4.48 eV. We assume that 
essentially all of this energy goes into rotational and vibrational excitation of the resulting 
H2 molecule, and hence is radiated away at low gas densities and is converted by collisional 
de-excitation into heat at high gas densities. 

Finally, we include hea ting due to the ionization of the gas by cosmic rays. Following 
Goldsmith fc Langerl (jl978l ). we assume that every ionization deposits 20 eV of heat in the 
gas, and so derive a heating rate that scales with the total cosmic ray ionization rate of 
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the gas (i.e. the sum of the ionization rates for the various individual chemical species, 
weighted by the fractional abundances of those species). Since considerable uncertainty 
exists concerning the value of the cosmic ray ionization rate in the local interstellar medium 
(see the discussion in iMcCall et al.ll2003l ). let alone concerning the appropriate rate to use 
in high-redshift protogalaxies, a more detailed treatment does not appear to be warranted 
at this time. 



4. Applications 

The chemical network and thermal model described in the preceding sections have a 
number of potential applications. One of the more obvious applications is the study of the 
cooling and gravitational collapse of gas in low metallicity protogalaxies. As previously noted, 
the density range for which our model is valid corresponds to a wide range of cosmological 
overdensities. It can therefore be used to study the thermal and chemical evolution of 
the majority of the gas within a given protogalaxy. For example, suppose we model the 
gas distribution within a z = 20 halo as a singular isothermal sphere with mean number 
density within the virial radius n = 0.4 cm~^ (a reasonable z eroth-order appro ximation for a 



protogalactic halo with a mean overdensity 6 = 200; see e.g. lAbel et al.ll2002l ). In this case, 
only ~4% of the gas within the halo has a density n > 100 cm~^ at which our model may 
break down. Therefore, although our model is of limited usefulness for studying gravitational 
fragmentation and star formation within this dense central region, it does allow one to study 
many other important problems, such as how the minimum protogalactic mass at which 
cooling becomes effective is affected by the presence of metals, or how the gas responds to 
the presence of an ultraviolet background. We use a version of this model to address some 
of these questions in paper II for the case of small protogalactic halos within initially ionized 
regions. 

A second possible application is the study of the evolution of the ISM in high redshift, 
metal-poor dwarf galaxies. In this case, much of the gas involved is often gravitational 
stable, and so the characteristic timescales are longer than the gravitational free-fall time. 
This limits the applicability of the model to lower densities than in the case of gravitationally 
collapsing protogalaxies. However, even if we adopt a relatively long characteristic timescale 
^char = 100 Myr, corresponding to the sound crossing time of a 1 kpc disk at a sound speed 
of 10 km s~^, the model remains valid for densities up to n ~ 5 cm~^ even in the least 
auspicious case (Z = 0.1 gas, with no significant UV radiation field present). Our model 
is therefore well-suited for use in the study of the evolution of the warm neutral component 
of the ISM in such galaxies and will in some cases also be useful in the study of the cold 
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neutral component. However, the reader is reminded that our model does not treat hot 
(T ^ 10"^ K), highly-ionized gas, and cannot be used to study the effects of stellar feedback 
(H II regions, superno yae, etc.) unless coupled w ith some existing model capable of treating 
this hot gas (see e.g. [Sutherland &: Dopit all 19931 ). 



As well as studying the chemical evolution of gas within dwarf galaxies and protogalaxies, 
we can also use it to study the chemical evolution of gas between galaxies, i.e. the intergalactic 
medium (IGM). One important reason to do so is the fact that several different obervational 
techniques have been s uggested that may allow one to probe the thermal and chemical state 
of this gas. lOhl (120021 ) has suggested that if large regions of the IGM remain neutral after 
the turn-on of the first observable ionizing sources, then atomic oxygen in the IGM may 
be detectable through its UV absorption, as i t will produce an O I forest analogous to the 
z < 6 HiLyman-a forest. iFurlanetto &: Loebl (120031 ) also argue that a substantial fraction 
of the metals in high-redshift Galactic winds are likely to be in low-ionization states and 
may be observable in absorption. Metals in the high redshift IGM may also leave detectable 



imprints in the fluctuation spectrum of the c osmic microwave background (|Basu et all 12004 



Hernandez- Monteagudo et al.l l2006l ) . Finally, iHernandez-Monteagudo et al.l (l2007l ) have re- 
cently suggested that fine structure emission from atomic oxygen may also produce a de- 
tectable spectral distortion of the CMB if the excited fine-structure levels can be populated 
by radiative pumping via the O I Balmer-a transition. Accurate modelling of many of these 
effects requires accurate modelling of the temperature and chemical make-up of the inter- 
galactic gas, and hence a model such as that presented here. 
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Table 1. List of the coUisional gas- phase reactions in our chemical model. 



No. 


Reaction 




Rate coefBcient (cm^ s ^) 




Ref. 


1 


H + e" ^ 


■H-+7 


fci = dex[-17.845 + 0.762 logT + 0.1523(logT)2 

- 0.03274(logT)3] 

= dex[- 16.420 + 0.1998(logT)2 

- 5.447 X 10"^(logT)'* 
+ 4.0415 X 10-5(logT)6] 


T < 6000 K 
T > 6000 K 


1 


2 


H- + H- 


H2 + e- 


k2 = 1.5 X 10-9 

= 4.0 X io-9T-o ir 


T < 300 K 
T > 300K 


2 


3 


H+H+ - 




ka = dex[-19.38 - 1.523 logT 
+ 1.118(logr)2 - 0.1269(logT)3] 




3 


4 


H + H+- 


H2 + H+ 


ki = 6.4 X 10-1° 




4 


5 


H- +H+ 


H + H 


fcs = 5.7 X 10-*^'r-o-5 + 6.3 X 10-^ 
- 9.2 X 10-"T0-5 + 4.4 X IQ-^T 




5 


6 


H++e-- 


H + H 


ke = 1.0 X 10-8 
= 1.32 X 


T < 617K 
T > 617K 


6 


7 


H2 + H+ 


-^H+ + H 


k7 = [-3.3232183 x IQ-^ 




7 



+ 3.3735382 x 10-^lnT 

- 1.4491368 X 10-^(lnT)2 
+ 3.4172805 X lO-^OnT)'^ 

- 4.7813720 X 10-^(\nT)^ 
+ 3.9731542 x 10-^"{lnTf 

- 1.8171411 X 10-"(lnr)6 



-39- 



Table 1 — Continued 



No. 


Reaction 




Rate coefficient (cm^ s ^) 




Ref. 


+ 3.5311932 X 10-i3(lnT)^] 

Xexp( -212^37. IS) 


8 


H2 + e- - 


H + H + e- 


ks = 3.73 X 10-9T0-1121 exp ( -%^^^) 




8 


9 


H2+H- 


> H + H + H 


fcg = 6.67 X 10-12TV2 exp [ (1 + ™)] 




9 


10 


H2 + H2 - 


H2 + H + H 


5.996 xlO'^^T*-^****^ „„„ f 54657.4 \ 
'^lO (1.0+6.761 xlO-<iT)5 ti««l exp J, ) 




10 


11 


H + e- ^ 


■ H+ + e" + e" 


fcii = exp[-3. 271396786 x 10^ 
+ 1.35365560 X lOMnTe 

- 5.73932875 x 10°(lnTe)2 
+ 1.56315498 x 10"(lnTc)^^ 

- 2.87705600 x 10"^(lnrc)'' 
+ 3.48255977 X lO-^llnT^f 

- 2.63197617 X 10-3(lnre)*^ 
+ 1.11954395 X lO-'^OnTe)^ 

- 2.03914985 x 10-6(lnre)S] 




11 


12 


D + e" ^ 


■ D+ + e" + e" 


fci2 = fell 






13 


H+ + e- - 




fcl3A = 1.269 X 10-13 (31'^614)1.503 

X [1.0+(604^625)0-470j_i.923 
fcl3B = 2.753 x 10-14 (31.^14)1-500 

X [1.0+ ("5^188 )0-4tf7] -2.242 


Case A 
Case B 


12 

12 


14 


D+ + e- ■ 


D + 7 


A;i4 = A;i3 







15 H- +e- ^ H + e- +e- fci5 = exp[-1.801849334 x 10^ 



11 
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Table 1 — Continued 



No. Reaction 



Rate coefficient (cm^ s ^) 



Ref. 



+ 2.36085220 x 10° 


InTe 


- 2.82744300 x 10' 


-1 


(lnTe)2 


+ 1.62331664 x 10" 


-2 


(lnTe)3 


- 3.36501203 x 10" 


-2 


(lnre)4 


+ 1.17832978 x 10" 


-2 


(lnTe)5 


- 1.65619470 X 10" 


^3 


(lnTc)6 


+ 1.06827520 x IQ- 


-4 


(InTe)^ 


- 2.63128581 x 10" 


-6 


(lnre)8] 



16 H- + H ^ H + H + e- kie = 2.5634 x iQ-^ry^^^^ 



Te < 0.1 eV 



11 



exp 
+ 1. 

- 1. 
+ 8. 

- 1. 
+ 2. 
+ 8. 
-2. 
+ 2. 



[-2.0372609 x 10^ 

13944933 x 10° InTe 



4210135 X IQ-^i 
4644554 x lO'^l 
4327641 X 10-3( 
0122503 X 10-4( 
6639632 x lO'^l 
5850097 X 10-5( 
4555012 X 10~S( 
.0683825 X lO'^l 



lnTe)2 
lnTe)3 

lnTe)4 

InTef 
InTe)^ 

lnTe)« 
lnTe)9] 



Te > 0.1 eV 



17 H- + H+ ^ + ( 



9T-0.35 



fci7 = 6.9 X lO-^T 



9.6 X 10-7T-0.90 



T < 8000 K 
T > 8000 K 



13 



18 H + D+ ^ D + H+ 



kis = 2.06 X 10-i°TO-396exp (- 
+ 2.03 X 10-9T-0-332 



33^ 
T ) 



14 



19 D + H+ ^ H + D+ 



A;i9 = 2.0 X 10-iOTO-402exp (-^) 
- 3.31 X 10-17T1-48 



T < 2 X 10^ K 



14 



= 3.44 X io-iot°-35 



T > 2 X 10^ K 
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Table 1 — Continued 



No. 


Reaction 




Rate coefficient (cm^ s ^) 




Ref. 




20 


H2 + D+ - 


HD + H+ 


k2o= [0.417 + 0.846 logT-0.137(logT)2] x lO'^ 




15 


21 


HD + H+ 


^ Ha + D+ 


fc2i = 1.1 X 10-9 exp(-^) 




15 


22 


H2 + D^ 


HD + H 


fc22 = 1.69 X 10-10 exp(-4f0) 
= 1.69 X 10-1° exp ( 46^80 + 


T < 200 K 

T > 200 K 


16 


23 


HD + H- 


*D + H2 


fc23 = 5.25 X 10-11 exp(-^) 

= 5.25 X 10-11 exp ( ^^o ^ 173900^ 


T < 200 K 
T > 200 K 


17 


24 


He + e" - 


^ He+ + e" + e" 


k2i = exp[-4.409864886 x lOi 
+ 2.391596563 X lOilnTe 

- 1.07532302 x lO^lnT^f 
+ 3.05803875 x lOOQnTo)^ 

- 5.6851189 X lO-i^nTe)" 
+ 6.79539123 x IQ-'^ilnT^f 

- 5.0090561 X 10-3(lnTe)S 
+ 2.06723616 x lO-i^nTe)^ 

- 3.64916141 X 10-6(lnTe)S] 




11 


25 


He+ + e- 


^ He + 7 


A:25,rr,A = IO-HT-"-"^ [12.72 - 1.615 logT 


Case A 


18 



- 0.3162(logT)2 + 0.0493(logr)3] 

fc25,rr,B = IQ-'^^T''^ '^ [11.19 - 1.676 logT Case B 18 

- 0.2852(logT)2 + 0.04433(logT)3] 

fc25 di = 1.9 X 10-3r-l-5 exp (-47^) 

X [l.0 + 0.3exp(-2ifM)] 19 
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Table 1 — Continued 



No. 


Reaction 




Rate coefficient (cm^ s ^) 




Ref. 




26 


He+ + H - 


He + H+ 


fee = 1.25 X 10-15 




20 


27 


He + H+ - 


He+ + H 


k27 - 1.26 X 10-9r-0-75 exp ( 12T500 ^ 
= 4.0 X 10-37^4.74 


T < 10000 K 
T > 10000 K 


21 


28 


He+ + D 


He + D+ 


fes = k26 




— 


29 


He + D+ 


He+ + D 


feg = fer 






30 


C+ + e- - 




feo = 4.67x10-12 (^^-0.6 

- 1.23 X 10-17 (g^j' ^'exp (21845.6^ 

= 9.62 X 10-8 exp (-115^786.2) 

\ 300 / ^ \ T / 


T < 7950 K 

7950 K < T < 21140 K 
T > 21140 K 


22 


31 


Si+ + e- • 


-^Si + 7 


fei=7.5xl0--(^)-°-^^ 
= 4.86 X 10-^2 -0.32 

= 9.08xl0--(^)°"' 


T < 2000 K 

2000 K < T < 10^ K 

T > 10* K 


23 


32 


0+ + e- - 


-^0 + 7 


fc32 = 1.30 X io-iot-0-64 

= 1.41 X iO'Wrj.-o-66 ^ 7 4 X -^Q-ij^-i.s 

xexp( ™^ [1 + 062 X exp ( ™)] 


T < 400 K 
T > 400 K 


24 


33 


C + e- ^ 


C+ + e- + e- 


fes = 6.85 X 10-8(0.193 + u)-iu0-25e-" 


u = 11.26/Te 


25 
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Table 1 — Continued 



No. 


Reaction 




Rate coefficient (cm^ s ^) 




Ref. 




34 


Si + e- ^ 


■ Si+ + e~ + e~ 


fc34 = 1.88 X 10-'^(1.0 + uO-5)(0.376 + uy^u'^-^^e''' 


u = 8.2/re 


25 


35 


+ e- ^ 


■ 0+ + e" + e" 


= 3.59 X 10-8(0.073 + u)-iu0-34e-« 


U = 13.6/Te 


25 


36 


0++H- 


^ + H+ 


he = 4.99 X io-iiTO-405 + 7.54 x iQ-Wrp-O-^^s 




26 


37 


+ H+ - 


^ 0+ + H 


ksr = [1.08 X 

+ 4.00 X 10-10T0 00669] exp {-^) 




27 


38 


+ He+ 


0+ + He 


fcss = 4.991 X 10 15 (^^^^J exp ( ^^^^^oo) 

■ r» Tror\ -ir\ — 15 / T \ —0.2163 / T \ 

+ 2.780 X 10 15 ( ^ ) exp ( ^J^^^ ) 




28 


39 


C + H+ - 


^C++H 


fc39 = 3.9 X 10-16T0-213 




27 


40 


C+ + H- 


^ C + H+ 


fc4o = 6.08 X 10-14 (,o^oj''%xp( ™) 




27 


41 


C + He+ - 


C+ + He 


kii = 8.58 X 10-17^0-757 

= 3.25 X 10-17T0-968 
= 2.77 X io-i9ri-597 


T < 200 K 

200 < T < 2000 K 

T > 2000 K 


29 


42 


Si + H+ - 


Si+ + H 


fc42 = 5.88 X 10-13TO-848 


T < 10* K 


30 



= 1.45 X lO-i^T T > 10* K 



-44- 



Table 1 — Continued 



No. Reaction 



Rate coefficient (cm^^ s ^) 



Ref. 



43 Si + He+ ^ Si+ + He = 3.3 x 10 



31 



44 C+ + Si ^ C + Si+ fc44 ^ 2.1 x 10" 



31 



45 Si+ + H+ ^ Si++ + H fc45 = 4.10 X 10-1° (j^) 



0.24 



1.0 + 3.17exp (^^3^)] exp (-^J^) 



30 



46 Si++ + H Si+ + H+ = 1.23 x 10-^ (j^)"-^"* 

X i.0 + 3.17exp(^^^) 



30 



47 Si++ + c- 



Si++7 



M7.rr 



= 1.75 X 10-12 (_T_) 



-0.6346 



fc47,di = 2.2552 X 10-"T-i-^exp (-^) 
+ 5.6058 X lO-^T-i-s exp (^-^) 



32 
33 



References. 



4: 
7: 



Karpas. Anicich &: Huntress ( 1979). 5: [Moseley et al 



I Savin et al 



Martin. Keogh fc Mandvl (Il998f). 11: 



Pou 



17: 



aert et al 



Wisliard (119791) 2: iLaunav et all (199lh. 3: iRamaker fc Peek 



200C 
199S 



24: 

28: 



( 1970l). 6: ISchneider et al 



(l2004h. 8 : IStibbe fc TennvsonI (Il999ir 9 : iMac Low fc ShuU 



Janev et al 



Zygelman et all (Il989h. 2 1: 



Nahar 



Zhao et al. 



Shavitt] Jl959l). 18 : iHummer fc StorevI (Il998r ). 19: rAldrovandi fc Peauignqtl Jl 973h. 20: 



"]l987l) 



12: 



119781), 14: ISavinI 15: I Cerlichl (Il982h. 16 



Ferland et al 



1976), 



19941 ). 



1986 ) 



Mielke et al 



Nahai 



iKimura et al (119931) . 22 : iNahar fc PradhanI (|1997), 23 
]l999l). 25 : IVoronovl Jl997l). 26: IStancil et all Jl999l). 27 : ^Standi et al 



( 20041) 29: IKimura et all (Il993f). 30 : iKingdon fc Ferland (1996|), 31: 



10: 
13: 



iW94D, 



Le Teuff. Millar fc Markwickl |2000h . 32: iNaharl |l995lll996l ). 33: iMazzotta et al.l (jl99 



Note. — T and Te are the gas temperature in units of K and eV respectively. References 
are to the primary source of data for each reaction. 
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Table 2. List of the photochemical gas-phase reactions in our chemical model. 



No. Reaction 



Cross-section (cm^) 



Reference 



H + 7^H++e- 



49 D + 7 ^ D+ 



50 He + 7 ^ Hc+ + e" 



CT48 = 6.3 X 10^18 (^) exp(4 - 4e-i arctan e) 
X [l-exp(-27r/e)]"^ 



0-49 _ cr4g 



-16 /'Bth\V2 



51 H-+7^H + c- 



ago = 3.1451 X 10-^" (^) 

1.0 - 4.7416 {^y^^ + 14.82 (-^) 

- 30.8678 {^f^^ + 37.3584 {^f 

- 23.4585 {^f^ + 5.9133 (p^f 

CT51 = 2.11 X io-if5(£; - £;th)^/2£;"^ 



52 H+ + 7 ^ H + H+ CT52 = dex 



-40.97+ 15.9795 f 

3 



3.53934 ( ^ 



+ 0.2581155 ( #■ 



dex 

+ 6.5785 X 10-2 



30.26 + 7.3935 f-^) - 1.29214 



53 H2 + 7 ^ H + H See ggXT] 

54 H2+7^H++e- (754 = 9.560 x lO^^^ - 9.4 x 10" 



17 



2.16 X 10-" ( #- 



1.51 X 10-1^ ( #- 



1.48 X 10-1^ 



-2.71 



55 HD + 7 ^ H + D 

56 C + 7 ^ C+ + e- 



See 

age = 5.027 x \^-^''F{x,y,y^,ya,P) 



Eth = 13.6 eV 

e = 



13.6 



1 



Eth 



13.6 eV 
24.6 eV 



Eth = 0.755 eV 

Eth = 2.65 eV 

2.65 <E < 11.27 eV 

11.27 <E < 21.0eV 



Eth = l5AeY 6 

15.4 <E < 16.5 eV 

16.5 <E < 17.7 eV 
17.7 < S < 30.0 eV 



Eth = 11.26 eV 
^ 1.133 



2.144 



X 

y = Vx^ + I.6O72 
2/^ = 0.09157 
2/a = 62.16 
P = 5.101 



57 + 7->0++e- a57 = 1.745 X 10-i5F(a;,y,y^,y„P) 



Eth = 13.62 eV 



Table 2 — Continued 



No. Reaction 



Cross-section (cm^) 



Reference 



1.240 



B.698 



y = Va:^ + 0.12712 
yy, = 0.07589 
ya = 3.784 
P = 17.64 



58 Si + 7 ^ Si^ 



ass = 2.506 x lQ-^^F{x, y, y^,ya, P) -Bth = 8.152 eV 



23.17 



1.672 X 10 



-5 



Va?2 + (I42072 



y 

y^ = 0.2837 
ya = 20.57 
P = 3.546 



59 



Si+ + 7 ^ Si^ 



"■59 



4.140 X 10-i8i;^(a;,y,2/^,j/„,P) 



16.35 eV 
^ 6.634 



2.556 



X — 

y = Va;2 + 0.12722 

= 1.570 
ya = 13.37 
P = 11.91 



References. 



Shapiro fc Kand (119871) 4: 



Osterbrock 



Wilms. Allen k McCravl (|2002[ ). 7: 



Dunn 



(1989) 
(1966), 



2: 
5: 



Yan. Sadeghpour fc Dalg arnol 1 1998h. 3: de Jone 
Draine & Bertoldi (1996), 6: lO'Neil fc Reinhardi 



Abgrall fc Romil (l2006h . 8: Iverner et all (|l99f}) 



(1972); 
|l978h : 



Note. — References are to the primary source of data for each reaction. E is the photon energy in eV and 
Eth is t he energy threshold in eV. The fitting function F used in the tabulated cross-sections for reactions 56-59 
is from lVerner et all |l996l ) and is given by P = [{x — 1)^ -|- yfu]y^'^^~^'^{\- + \fyJvo)^^ ■ Photodissociation of 
H2 and HD occurs via absorption into a large number of discrete spectral lines and so no simple cross-section 
can be given for these processes; see ii2.4.1l for more details 
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Table 3. List of the grain surface reactions included in our chemical model. 



No. Reaction 



Rate coefficient (cm s 



3 



Ref. 



60 H + H ^ H2 

61 H+ + e- 

62 D+ + e" 

63 He+ + e" 

64 C+ + e- 

65 0+ + e- 

66 Si+ + e- 



C 

O 
Si 



ho = 3.0 X 10-^^T^-^{V/Vq)[1.0 + 4 X 10-2(T + T^^f-^ 



D k. 



62 



65 



1.0 + lO^exp - 



+ 2 X 10-^T + 8 X lO^'T 



H kai = 1.225 x 10-^''^{V/Vq)[1.0 + 8.074 x lO^^^^-^^^ 

(1.0 + 5.087 X 102r0-01586^-0.4723-1.102xl0-Mnr)]-l 



V2 



He fees = 5.572 x 10-^^{V/Vq) [1.0 + 3.185 x 10"^^ 



7„/,1.512 



:i.0 + 5.115 X lO^T^' 



903x10^'^ ^/,-0. 4956-5. 494x 10~^ In T\l -1 



= 4.558 x 10-i3(p/2)q)[1.0 + 6.089 x lO-^^^-^^s 
(1 + 4 331 X io^T°'°^^'*^-?/'~°''^^^°~"'^'^'^^^^°~*^^-'")]~^ 



1^61 



kee = 2.166 x lO-i^(r'/r'0) [1.0 + 5.678 x 10"^ 



1.874 



'1 + 4 375 X lQ^T^-^'^^x'^0^^^-0-89M-7. 538x10--' InT-s^j-l 



3 
2 



3 

2 



Note. — P is the dust-to-gas ratio and Vq is the dust-to-gas ratio in the local ISM. We gen- 
erally assume that V/Vq = Z/Zq. T and Tgr are the gas and grain temperatures, respectively. 
The parameter ip in the grain recombination rates is given by = CVT/n^, where G - O.OlJ^i 
i s a rn easure of the radiation energy density between 6 eV and 13.6 eV relative to the iHabing 
f|l968h field. 



References. — 1: iHoUenbach fc McKed fll979r ): 2: IWeingartner fc Draind fl2001af ): 3: This 
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work, but based on IWeingartner fc Draind (l2001al ) 



Table 4. List of cosmic ray ionization processes included in our chemical model 



No. Reaction 



Ci/Cn 



Ref. 



67 H + c.r. ^H+ + e- 1.0 Le Teuff. Millar fc Markwick r200nl 

68 D + c.r. ^ D+ + e" 1.0 Le Teuff. Millar fc Markwick QOOO) 

69 He + c.r. ^ He+ + e" 1.09 Le Teuff. Millar fc Markwick r2000) 

70 B.2 + c.r. H+ + e" 2.0 Le Teuff. Millar fc Markwick (20001 

71 C + c.r. ^C+ + e- 3.83 Le Teuff. Millar fc Markwick (20001 

72 O + c.r. 0+ + e" 5.67 Le Teuff. Millar fc Markwick (20001 

73 Si + c.r. Si+ + e" 6.5 Lotz (19671 : Langer (1978 1 

74 Si+ + c.r. ^ Si++ + e- 2.5 Lotz (19671 : Langer (19781 



Note. — We list here the ratio of the various rates to the rate of process 
67, the cosmic ray ionization of atomic hydrogen, Ch? which we treat as an 
adjustable parameter in our models. Rates for co smic ray ioniz ation of Si 
and Si+ we re cal c ulated following the prescription in iLangen (Il978l 1 and using 
data from iLotzl (1l967l 1 under the assumption that the effective number of 
outer shell electrons for Si and Si"*" in the high energy limit is the same as 
that for C and C+. 
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Table 5. Atomic data for the fine structure transitions included in our thermal model 

Coolant Transition gj Ajj(/xm) Eji/k (K) Aj^ (s~^) 

C 1^ 3 1 609.2 24 7.9 x 10"^ 

C 2^ 5 1 229.9 63 2.1 x 10"^^ 

C 2^ 1 5 3 369.0 39 2.7 x 10"^ 

O 1^ 3 5 63.1 230 8.9 x 10"^ 

O 2^ 1 5 44.2 330 1.3 x 10"^° 

O 2^ 1 1 3 145.6 98 1.8 x 10"^ 

Si 1^ 3 1 129.6 110 8.4 x lO"*^ 

Si 2^ 5 1 44.8 320 2.4 x 10"^° 

Si 2^ 1 5 3 68.4 210 4.2 x 10"^ 

C+ 1^ 4 2 157.7 92 2.3 x lO^^ 

Si+ 1^ 4 2 34.8 410 2.2 x 10"^ 
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Table 6. CoUisional de-excitation rates for atomic fine-structure coolants 



Coolant Collider 



De-excitation rates (cm^ s ) 



Temperature range (K) Refs. 



C 



o-H, 



C 



P-H2 



gio = 8.7 X 10-" - 6.6 x 10"" exp {-^^ 

+ 6.6 X 10-"exp(-2fo) 
(720 = 1.2 X 10-1° - 6.1 X 10-" exp 
2.9 X 10-1° - 1.9 X 10-1° 



921 



exp 



910 = 7.9 X 10-" - 8.7 X 10-" exp 
+ 1.3 X 10-10 exp (-^) 

q2o = 1.1 X 10-1° - 8.6 X 10-11 exp 
+ 8.7X 10-11 exp (-Jl^) 

921 = 2.7 X 10-1" - 2.6 X 10-1° exp 
+ 1.8 X 10-1° exp (-j;^) 



387.3/ 

I—) 

348.9/ 

126.4 y 

Z_) 

223.0/ 

I—) 

250.7/ 



c 


H 


910 = 1.6 X 10-i°T2° i^ 

920 = 9.2 X 10-iir2°-26 

921 = 2.9 X 10-i°T°-26 




2 

2 
2 


c 


H+ 


910 = (9.6 X 10-11 - 1.8 X 10-"T-|- 1.9 x io-i»T^)T'^ 'i5 
= 8.9 X io-iot° ii7 


T < 5000 
T > 5000 


3 
3 






920 = (3.1 X 10-12 - 6.0 X 10-i^T -|- 3.9 x 10-^°T^)T 
= 2.3 X io-9t° °965 


T < 5000 
T > 5000 


3 
3 






921 = (1.0 X 10-1° - 2.2 X lO-i-iT + 1.7 x io-i&T^)To ro 
= 9.2 x io-9t° °535 


T < 5000 
T > 5000 


3 
3 


c 


e- 


910 = 2.88 X 10-6t-°-5 exp[-9. 25141 - 7.73782 x lO-i InT 
+ 3.61184 X 10-i(lnr)2 - 1.50892 x 10-2(lnr)3 
-6.56325 x 10-^(lnT)4] 


T < 1000 


4 






= 2.88 X 10-6T-°-5exp[-4.44600 x 10^ - 2.27913 x lO^lnT 
+ 4.2595 x 10i(lnr)2 - 3.47620 x 10°(lnr)3 
-M.0508 X 10-i(lnT)4] 


T > 1000 


4 



920 = 1.73 X 10-6T-°-5exp[-7.69735- 1.30743 InT 
+ 0.697638(lnr)2 - 0.111338(lnr)3 
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Table 6 — Continued 



C()()l;uit C\)lli(loi; D(^-oxcit;itiou ralos (cm'' s ' ) T(nu])(>ratur(' raiigi." (K) Rcfe. 



+ 0.705277 X 10-2(lnT)4] T < 1000 4 

= 1.73 X 10"6j.-o.5^j.p[3 50609 x 10^ - 1.87474 x lOMnT 
+ 3.61803 X 10i(lnr)2 - 3.03283 x 10°(lnr)3 

+ 9.38138 x 10-2(liiT)4] T > 1000 4 

921 = 1-73 X 10-6^-0.5 gj.p[_7 4387 _ 0.57443 InT 
+ 0.358264(lnT)2 - 4.18166 x 10-^{\nTf 

+ 2.35272 x lO-^'OnT)*] T < 1000 4 

= 1.73 X IQ-^T-^-^ exp[3.86186 x 10^ - 2.02192 x 10^ InT 
+ 3.85049 x 10i(lnT)2- 3.19268 x 10°(lnT)3 







+ 9.78573 X 10-2(lnT)4] 


T > 1000 


4 





0-H2 


gio = 2.7 X io-iiTO-362 




5 






920 = 5-49 X l0-iir" -'^i7 




5 






921 = 2.74 X 10-14^1060 




5 





P-H2 


910 = 3.46 X io-iiTO-316 




5 






920 = 7.07 X io-iiTO-268 




5 






921 = 3.33 X io-i5ri-360 




5 





H 


910 = 9.2 X 




5 






920 = 4.3 X lO-iiTO-80 




5 






921 = 1.1 X IO-iOtO-^^ 




5 





H+ 


910 = 6.38 X lO-iiTO -io 


T < 194 


6 






= 7.75 X 10-12TO-80 


194 < T < 3686 








= 2.65 X lo-iorO-37 


T > 3686 








920 = 6.10 X io-i3t1i° 


T < 511 


6 






= 2.12 X io-i2r0-90 


511< T < 7510 








= 4.49 X io-iotO-30 


T > 7510 








921 = 2.03 X io-iir"-56 


T < 2090 


6 






= 3.43 X 10-10T019 


T > 2090 
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Table 6 — Continued 



Coolant 


Collider 


De-excitation rates (cm^ s 




Temperature range (K) 


Refs. 





e" 


(710 = 5.12 X I0-i"r-o o75 
(720 = 4.86 X io-i"T-0-026 

921 = 1.08 X 10-14^0-926 






7 
7 
7 


Si 


H 


910 = 3.5 X 10-ioT2-° °^ 

920 = 1.7 X 10-11T017 

921 = 5.0 X lO-iOTsO-i^ 






2 
2 
2 


Si 


H+ 


9io = 7.2 X 10-9 

920 = 7.2 X 10-9 

921 = 2.2 X 10-8 






2 
2 
2 


C+ 


0-H2 


910 = 4.7 X 10-10 + 4.6 x 10" 
= 5.85 X 10-10T0 07 


-13ji 


T < 250 
T > 250 


8, 9 


c+ 


P-H2 


910 = 2.5 X lo-iOTO-12 
= 4.85 X io-iotO-07 




T < 250 
T > 250 


8, 9 


c+ 


H 


910 = 8.0 X lO-iOTsO-O^ 
= 3.1 X lO-iOrO-385 




T < 2000 
T > 2000 


2, 10 


c+ 


e" 


910 = 3.86 X 10-^T2-°-^ 
= 2.43 X 10-^T2-°-^''^ 




T < 2000 
T > 2000 


11 


Si+ 


H 


910 = 4.95 X 10-iOr20-24 






12 


Si+ 


e~ 


910 = 1.2 X 10-^X2°-^ 






13 



Note. — 



0-H2 and P-H2 denote ortho-H2 and para-H2 respectively. T is the gas temperature 
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(in Kelvin) and T2 = IQ-^T. 



References. 



1: ISchroder et aD (Il99lh: 2: iHollenbach fc McKeel (Il989l) : 3: 



Roueff fc Le BourlotI (llQQOl ) : 4: I Johnson et alj (119871): 5 : Flow er, private communica t ion; 6 : 
PeauignotI |l99d Il996l )~ bell. Berrington fc Thonia] |l998h : 8: iFlower fc Launavl |l977l ): 



9: assumed to have t he same scaling with T a s the l ow temperatur e H ra te fo r temperatures 
above the range of the'Fl ower fc LaunavlJl977l) fit: 10: IXeenan et al\ |l986h : 11: Iwilson fc Bell 
I2OO2I ): 12: ,Roueft (,199Q) : 13: lOufton fc KingstonI ilMll), extrapolated to T < 4000 K assum- 
ing constant collision strength 
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Table 7. Other processes included in our thermal model. 

Process Rate (erg cm~^ s~^) Ref. 

Cooling: 

H excitation A = 7.5 x IQ-^^ (^1.0 + a/T/IQs) \xp {-^^^) n^nn 1 

He excitation (l^S state) A = 1.1 x 10-^^T^-0S2 g^p (- ™0) ^eUHe 2 

He excitation (23s state) A = 9.1 x 10-27^-0.1687 J^^ g + ^/T/W^^ \xp {-^^) nln^^+ 1 
He+ excitation A = 5.54 x 10-17^-0.397 q _^ y/T/W^ ^ exp (-^^f^) nen^^+ 1 

H collisional ionization A = 2.179 x W^^^ knUenn 3 

He collisional ionization A = 3.94 x 10~^^k24nen}ie 3 

Compton cooling A = 1.017 x IQ-^'^T^^^ (T - Tcmb) 1 

Bremsstrahlung A = 1.426 x IQ-^"^ ZfT^^^gs{Zi,T)neni 4 

gs = 0.79464 + 0.1243 log (T/Zf) {T/Zf) < 320000 K 

= 2.13164 - 0.1240 log {T/Zf) [T/Zf) > 320000 K 

H^ recombination (radiative) A = 1.38 x 10^^^TA;i3nenH+ 5 

He"*" recombination (radiative) A = 1.38 x 10~i^TA;25,rrf^e'T'He+ ^ 

He^ recombination (dielectronic) A = 6.54 x 10^^"'^A;25,dif^e'^He+ 

Grain surface recombination A = 2.33 x iQ-^Orp0.9i^o.7i/T'^ °'^^ (^) ^ 

H2 rovibrational lines See ^3.21 9 

HD rovibrational lines See ^3.21 10 
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Table 7 — Continued 



Process Rate (erg cm s ^) Ref. 

H2 collisional dissociation A = 7.2 x lO^"*^^ (^g^i-H + kiQU-n^) 11 

Gas-grain energy transfer A = 3.8 x 10-33^1/2^^ - Tg,) [l.O - O.Sexp (-f )] 12 
Heating: 

Pliotoelectric effect T = 1.3 x W'^^e G (^^) n 13 

^ 4.9x10-2 , 3.7xlO-2(T/10000)0-7 

1.0+4.0x10-3^0.73 "T 1.0+2.0x10-4?/; 

H2 photodissociation F = 6.4 x l^^^^ Rd\ssn-B_2 14 

UV pumping of H2 T = 2.7 x lO-iii^dissnn^ (^) 15 

H photoionization Dependent on incident spectrum; see ^2.41 16 

He photoionization Dependent on incident spectrum; see ^2.41 17 

Gas-phase H2 formation T = 2.93 x lO-^^/j^nn- + 5.65 x lO^^^/.^^^^ {l^T^^ 

H2 formation on dust grains F = 7.16 x 10~^^/i;6o renn ( ^ i!l ) 19 

Cosmic-ray ionization F = 3.2 x lO^^^Ctot'^' 20 



References. — 1: ( 


Der 


1 (1992). 2: 


Brav et al. (2000) 


3: 


Janev et al 


• ( 


1987 


). 4: Shaoiro k, Kane 


(1987) 5: Ferland et al. (19921). 6: Hummer & Storev ( 


1998 


). 7: 


Aldrovandi & Peauienot ( 


1973). 


8: 


Wolfire et al 


(I2OO3). 9: 


Le Bourlot et al. ( 


1999). 10: LiDovka. Niihez-LoDez. & Avila-Rees( 
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2005 


), 
), 


11: 
13: 


Mac Low & Shul 


(1986): Martin. Keoeh & Mand^ 


^ (Il998). 12: ] 


Hollenbach k. McKee ( 


1989 


Bakes &: Tielens 


I994I): Wolfire et al. 


(1995). 14: 


Black &; Dalgarnc 


(1977). 15: Burton et al. 


(199C). 16: 


Osterbroc 


k ( 


19891). 17: " 


Y^an. Sadeghpour &; Dalgarnc 


(1998). 


18: 


Launav et al. ( 


1991 
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Karpas. Anicich fc HuntressI (|l979l ) 
(Il97i 



19: iHollenbach k McKed (|l979l ^. 20: iGoldsmith k Langei 



Note. — Zi and Ui are the ion charge and number density of ion i. The parameter ip is given by 
tjj = Gy/T/O.brLe, whe re G ~ O .OlJo i is a measure of the radiation energy density between 6 eV and 
13.6 eV relative to the lHabingI (119681 ) field. .Rdiss is the photodissociation rate, calculated as discussed 
in §2.4. 1[ ^tot is the total cosmic-ray ionization rate (i.e. the sum of the rates for the various different 
species, weighted by their fractional abundances: Ctot = Yli ^id)- Finally, note that our treatment of 
recombination cooling here is approximate, but that it should be accurate enough for most purposes. 



